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SUMMARY 


In  this  report,  the  state  of  the  technology  in  shell  nosing  and  the 
theoretical  and  experimental  studies  on  shell-nosing  process  are  briefly  reviewed. 
In  order  to  determine  an  optimum  combination  of  the  process  variables,  a  number  of 
computerized  mathematical  models  for  cold  and  hot  nosing  of  shells  were  developed. 
All  these  models  were  assembled  in  one  comprehensive  computer  programm  called 
NOSING. 

Various  capabilities  for  computerized  modeling  of  the  nosing  process, 
developed  in  this  project,  are: 

•  Preform  Design.  Based  on  the  consideration  of  local  strains 
in  the  deformed  shell  during  nosing,  a  method  for  designing 
of  nosing  preforms  was  computerized.  This  method  for  preform 
design  is  capable  of  considering  elongation  due  to  nosing,  and 
it  has  also  been  programmed  on  HP-67  programmable  calculator. 

•  Temperature  Distribution  Prior  to  Nosing.  A  mathematical  model 
to  predict  time-dependent  temperature  distributions  due  to 
induction  heating  of  preform  prior  to  hot  nosing  was  developed 
and  computerized.  The  analysis  assumes  uniform  heat  generation 
along  the  length  of  the  tube  inside  the  coil  and  neglects  the 
end  effects.  Therefore,  the  predictions  are  less  accurate  near 
the  end  of  the  induction  coil. 

•  Prediction  of  Metal  Flow  During  Nosing.  A  mathematical  model 
for  predicting  metal  flow  in  nosing  of  shells  was  developed. 

This  model  considers  preforms  with  uniform  wall  thickness. 

Thus,  the  model  is  exactly  valid  for  nosing  of  shells  up  to 
105-mm  size  shells,  which  are  normally  nosed  from  preforms 
with  uniform  or  near-uniform  wall  thicknesses.  For  larger 
shells,  where  the  preform  wall  thickness  is  not  uniform,  this 
analysis  can  be  applied  only  approximately. 

•  Load-Stroke  Curve  in  Nosing.  In  order  to  generate  the  load- 
stroke  curve  in  the  nosing  operation,  using  preforms  with 
non-uniform  wall  thickness,  a  computer  program  to  simulate 
the  nosing  process  was  developed.  This  program  simulates  the 
nosing  process  in  a  finite  number  of  discrete  steps  and 
utilizes  Nadai's  stress  analysis,  and  considers  the  flow  stress 
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of  the  deforming  material  as  a  function  of  the  strain, 
strain  rate  and  temperature.  It  also  checks  for  local 
bulging  at  the  nose  base,  or  for  Euler's  buckling,  at 
each  step  of  simulation. 

•  Computer  Program  NOSING.  Finally,  all  the  mathematical 
models  described  above  were  integrated  in  a  comprehensive 
computer  program  named  NOSING.  This  program  is  capable 
of  simulating  both  cold  and  hot  nosing  of  shells  and 
determining  the  optimum  combination  of  process  parameters. 


INTRODUCTION 

In  all  modern  methods  of  shell  manufacturing,  the  cavity  is  formed  to 
finish  shape  and  the  machining  is  restricted  to  the  outer  surface  of  the  shell. 

The  open  end  of  the  rough-machined  shell  is  closed  in  and  the  ogive  is  formed. 

The  closing  in  is  accomplished  by  forcing  a  contoured  die  axially  over  the  open 
end  of  the  shell,  while  the  body  of  the  shell  is  well  supported  by  a  chuck,  as 
shown  in  Figure  1.  In  all  high-explosive  shells,  from  75  mm  to  240  mm,  the  fuze 
thread  diameter  is  the  same.  Therefore,  the  open  end  of  the  largest  shell  must 
be  deformed  about  three  times  as  much  as  the  smaller  calibers  to  produce  the 
same  size  fuze  hole.  As  a  result,  155-mm  shells  and  larger  sizes  are  hot  nosed, 
while  the  75-mm  to  105-mm  shells  can  be  cold  nosed. 

The  present  program  was  aimed  at  increasing  the  productivity  in  nosing 
of  artillery  shells.  The  flow  of  metal  in  nosing  is  very  complex  and  a  slight 
variation  in  the  friction  or  the  temperature  conditions  may  result  in  misformed 
shell  due  to  improper  metal  flow.  There  is  very  little  quantitative  technical 
information  available  and,  in  most  shell-manufacturing  plants,  extensive  experience 
and  expensive  trial-and-error  techniques  are  necessary  in  order  to  design  the 
nosing  process.  Further,  cold  or  hot  nosing  of  shells  is  traditionally  carried  out 
in  hydraulic  presses  at  moderate  speeds.  Recently,  high-speed  mechanical  presses 


J 


4 


are  being  considered  for  the  nosing  operation.  The  process  design,  as  practiced 
today,  leaves  much  to  experience  and  intuition.  In  such  situations,  mathematical 
modeling  of  the  process  capable  of  predicting,  for  example,  conditions  leading  to 
buckling  in  hot  nosing  or  extent  of  wall  thickening  in  cold  nosing  are  essential 
in  order  to  eliminate  expensive  trial-and-error  methods.  Thus,  an  objective  and 
reliable  procedure  can  be  developed  to  select  the  optimum  equipment,  tooling,  and 
operating  conditions. 


OBJECTIVES 


The  overall  objective  of  the  present  program  was  to  develop  analytical 
techniques  to  determine  the  optimum  combination  of  process  variables  for  defect- 
free  nosing  of  shells.  Both  cold  and  hot-nosing  operations  were  considered. 

The  specific  objectives  of  the  present  program  were  the  following: 

(a)  In  hot  nosing,  determine  temperature  distribution  in 
the  shell  wall  due  to  preheating. 

(b)  Develop  a  method  for  predicting  the  load-stroke  curve 
in  both  cold  and  hot  nosing  under  a  specified  combina¬ 
tion  of  process  variables,  such  as  die  configuration, 
wall  temperature,  lubricant,  and  speed  of  operation. 

(c)  Develop  a  criterion  for  predicting  buckling  due  to 
axial  loads  in  hot  and  cold  nosing. 

(d)  Determine  the  effect  of  ram  speed  on  the  forming  load, 
metal  flow,  and  possibility  of  buckling  in  both  cold 
and  hot  nosing. 

(e)  Using  the  information  developed  in  the  above  items  (a) 
through  (d) ,  develop  a  computer  program  to  determine 
optimum  combination  of  process  variables  for  cold  and 
hot  nosing  of  shells. 
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PROGRAM  HIGHLIGHTS 


The  present  program  was  conducted  over  a  period  of  18  months  and  the 
work  under  it  included  the  following  major  tasks: 

Task  1.  Temperature  Distribution  after  Preheat 

The  coils  in  the  preheaters  provide  the  heat  source  to  the  shell  during 
nosing.  By  modeling  the  coils  as  an  energy  source,  the  heat  conduction  equations 
for  the  shell  were  solved.  The  solution  technique  provided  the  temperature 
distribution  in  the  shell  as  a  function  of  location  and  time.  Cool-off  due  to  any 
delays  between  preheating  and  actual  nosing  was  also  accounted  for. 

Task  2.  Metal  Flow  and  Velocity  During  Nosing 

The  distribution  and  velocity  of  the  metal  flow  was  modeled  using  one 
of  the  methods  of  plasticity.  Temperatures  (from  Task  1),  stroke,  strain  rate, 
friction  and  die  geometry  were  included  in  the  metal-flow  calculations. 

Task  3.  Stress  Analysis 

A  numerical  extension  of  the  original  Nadai  stress  analysis  in  nosing 
was  computerized.  The  results  were  an  estimate  of  the  axial,  radial  and  circum¬ 
ferential  stresses  and  strains  during  nosing. 

Task  4.  Material  Failure  Analysis 

Using  the  stress  analysis  (Task  3),  a  criterion  for  local  plastic 
buckling  or  bulging  was  formulated.  This  criterion  contains  the  effects  of 
friction,  temperature,  strain  and  strain  rate. 
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Task  5.  Optimization  Computer  Program 

A  computer  program  utilizing  the  results  of  Tasks  1  to  4  was  written. 
Parametric  study  to  determine  optimal  nosing  conditions  of  selected  shells  was 
performed. 


BACKGROUND  ON  SHELL  NOSING 


The  small  caliber  shells  (up  to  105  mm)  can  be  cold  nosed  while  larger 
shells  (155  mm  and  up)  are  hot  nosed.  In  cold  nosing,  initially  the  whole  shell 
is  at  constant  temperature  and  it  has  uniform  flow  stress.  However,  the  deforma¬ 
tion  work  hardens  the  shell  material  nonuniformly ,  and  the  resistance  to  further 
deformation  by  various  parts  of  the  shell  is  also  unequal.  The  elastic  spring 
back  is  very  small  compared  to  plastic  strain  and,  therefore,  the  dies  are  given 
the  shape  of  rough-finished  nose.  The  slight  elastic  back  spring  is  removed  in 
the  process  of  finish  turning. 

During  the  nosing  operation,  the  shell  metal  in  the  die  can  move  either 
radially  inward,  or  in'  the  direction  of  tangent  to  the  die  profile,  or  in  both  the 
directions.  Since  metals  under  plastic  state  flow  in  the  direction  of  least 
resistance,  the  flow  of  metal  in  cold  nosing  is  affected  by  the  characteristics 
of  the  lubricant  and  the  speed  of  the  nosing  operation.  If  the  lubrication  is  not 
sufficient  and  if  the  nosing  is  slow,  the  lubricant  is  squeezed  out  and  the  friction 
between  shell  and  die  not  only  keeps  the  shell  from  lengthening,  but  shortens  it. 
Conversely,  if  the  lubrication  is  adequate  and  if  the  nosing  is  done  rapidly,  the 
lubricant  is  squeezed  out  to  a  lesser  extent;  the  friction  is  reduced  and  the  shell 
is  lengthened.  The  finish  of  the  external  shell  surface  in  rough  turning  also 
affects  the  metal  flow  in  nosing;  smoother  finish  encourages  the  escape  of  the 
lubricant  while  very  rough  surface  finish  results  in  compression  of  the  ridges 
left  from  machining  and  provides  improved  lubrication. 

In  hot  nosing,  the  open  end  of  a  forged  shell  is  preheated  radially  to 
temperatures  between  1500  F  to  1900  F  by  tubular  wound  induction  coils.  Then  the 
shell  is  nosed  by  forcing  it  into  a  suitably  shaped  die.  The  flow  stress  of  the 
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deforming  material  under  these  conditions  is  a  function  of  strain  rate  and 
temperature.  The  temperatures  in  the  shell  wall  are  influenced  by  (a)  preheating 
conditions,  (b)  heat  generation,  due  to  deformation  and  friction,  and  (c)  heat 
transfer  to  the  nosing  die.  Determination  of  these  temperature  distributions 
along  the  length  of  the  shell  wall  is  important  and  necessary  in  order  to  predict 
accurately  the  axial  nosing  loads. 

Because  of  the  temperature  distribution  along  the  axis  of  the  shell, 
the  flow  of  metal  in  hot  nosing  is  also  complex.  In  fact,  the  temperature  distri¬ 
bution  is  rather  critical  and  much  experimentation  has  been  done  on  the  heating  and 
on  the  contour  of  rough-turned  shells.  Preheating  is  done  at  a  certain  rapid  rate 
so  that  the  high  temperature  does  not  travel  back  to  that  part  of  the  shell  which 
should  remain  comparatively  cold;  otherwise,  buckling  of  the  shell  wall  occurs,  as 
shown  in  Figure  2. 


Fig  2  Buckling  at  the  Nose  Base  in  Hot  Nosing  of  Shells 


Lubrication  plays  an  important  role  in  hot  nosing  also.  Usually,  the 
lubricant  is  sprayed  on  rough-turn  shells  before  preheating  for  the  nosing 
operation.  The  lubricant  evaporates  somewhat  unevenly  during  preheating  and  a 
nonuniform  lubrication  may  exist  prior  to  nosing.  This  causes  an  uneven  or 
oblique  top.  Further,  lack  of  adequate  lubrication  causes  the  shell  to  be  short 
and  the  walls  to  become  thicker,  whereas  good  lubrication  results  in  longer  shells 
with  less  wall  thickening.  The  influence  of  die  temperature  and  the  speed  of 
nosing  further  complicates  the  situation. 

The  Variables  of  the  Shell-Nosing  Process 

As  discussed  earlier,  many  factors  influence  the  shell-nosing  process. 

In  order  to  optimize  the  process,  the  variables  which  determine  the  process 
conditions  must  be  considered  individually  and  collectively,  and  their  contribution 
to  the  process  must  be  defined  quantitatively.  The  significant  variables  of  the 
shell-nosing  process  are  listed  in  the  following: 

(A)  Projectile  Material  Variables 

•  Flow  stress  and  its  dependency  on  strain,  strain 
rate  and  temperature 

•  Physical  and  thermal  properties  (density,  specific 
heat,  thermal  conductivity) 

•  Metallurgical  properties 

(B)  Tooling  Variables 

•  Die  configuration 

•  Die  and  container  materials  and  their  properties 

(C)  Process  Variables 

•  Projectile  forging  dimensions 

•  Depth  of  penetration  in  the  die 

•  Projectile  forging  temperature 

•  Die  temperature 

•  Speed  of  operation 

•  Lubrication  and  the  lubricant  properties 

•  Nosing  load  and  the  die  pressure 


•  Limiting  load  and  penetration  for  buckling 

•  Variations  in  the  wall  thickness  and  length  of 
the  projectile 

•  Alignment  of  nosing  die  and  nest 
(D)  Product  Properties 

•  Concentricity  of  projectile,  free  from  bulging 
and  buckling 

•  Dimensional  tolerances  and  surface  finish 

•  Mechanical  properties  (tensile,  fatigue,  etc.) 

•  Metallurgical  properties  (grain  size,  phase 
transformation,  etc.) 


A  Brief  Review  of  Theoretical  and  Experimental 
Studies  on  Shell  Nosing 

Nadai  (Ref .  1)  conducted  an  extensive  theoretical  investigation  on  the  forces 
required  in  nosing  of  shells  for  the  ASME  Special  Research  Committee  on  Forging  of 
Steel  Shells.  Although  hot  nosing  was  of  primary  interest,  Nadai' s  developments 
can  be  extended  to  the  cold-nosing  process  too.  This  work  is  basically  an  extension 
of  the  theory  of  curved  shells  to  cover  cases  in  which  the  metal  is  in  a  plastic 
state  of  equilibrium.  The  theory  is  based  upon  constant  coefficient  of  friction  and 
ignores  variations  in  the  flow  stress  of  the  shell  material.  General  equations  for 
both  conical  and  curved  shell  nose  were  established  for  the  plastic  state  and  for 
variable  wall  thickness.  However,  investigation  with  uniform  wall  thickness  was 
considered.  The  distribution  of  the  meridianal  and  circumferential  stress  and  of  the 
nosing  pressure  were  studied  for  a  number  of  cases.  Nadai  also  attempted  to  analyze 
the  distortion  of  the  metal  elements  during  the  formation  of  the  nose,  and  develop 
expressions  to  predict  the  original  contour  of  the  shell  which  after  nosing  furnished 
a  prescribed  profile  on  the  nose  of  the  shell.  Onat  and  Prager  (Ref . 2)  extended  Nadai's 
work  and  Included  the  changes  in  the  shell-wall  thickness  due  to  the  nosing  operation. 
They  investigated  the  influence  of  these  changes  on  the  stresses  in  nosing.  A 
linearized  theory  of  nosing  of  shells  has  been  presented  by  Singh  (Ref. 3)  who  has  shown 
that  excellent  approximations  to  the  predictions  of  von  Mises'  theory  can  be 
obtained.  However,  none  of  these  analyses  take  into  account  (a)  the  dependence  of 
the  flow  stresses  on  strain  in  cold  nosing,  and  on  strain  rate  and  temperature  in 


hot  nosing,  (b)  variation  in  temperature  along  the  length  of  the  shell  nose  in 
hot  nosing,  and  (c)  variation  of  the  friction  coefficient  in  the  axial  direction. 

The  published  literature  on  nosing  of  shells  also  includes  a  few 
systematic  experimental  studies.  The  model  tests  by  Carlson  (Ref . 4) ,  conducted  for 
the  ASME  Special  Research  Committee  on  Forging  of  Steel  Shells,  show  the  effects 
of  temperature,  lubrication,  and  speed  on  the  nosing  process.  Carlson  has  also 
presented  an  analysis  of  the  strain  distributions  and  has  indicated  a  method  for 
predicting  the  nose  profile  necessary  to  produce  a  desired  finish  shape.  Recently, 
Cruden  and  Thomson  (Ref .  5)  conducted  an  experimental  study  of  the  nosing  process  to 
establish  the  limitations  of  the  process  and  to  assess  the  effects  of  the  various 
process  parameters.  Some  of  the  practical  aspects  of  shell  nosing,  based  on 
surveys  of  plant  practices,  are  summarized  in  a  report  by  Veth(Ref.6),  et  al.  However, 
the  recommendations  of  these  studies  are  largely  qualitative  in  nature  and  cannot 
be  used  reliably  for  process  optimization  without  an  extensive  investigation  of 
similar  nature. 


PREFORM  DESIGN  FOR  NOSING  OF  SHELLS 


During  the  nosing  operation,  the  shell  metal  inside  the  die  can  move 
either  radially  inward,  or  in  the  direction  of  the  tangent  to  the  die  profile,  or 
in  both  directions.  Certain  shell  specifications  require  uniform  wall  thickness 
after  nosing,  such  as  in  the  case  of  175  mm,  M437  shell;  whereas,  in  the  case  of 
other  shells,  such  as  155  mm  M107,  the  required  wall  thickness  after  nosing  is 
defined  by  two  ogive  radii,  as  shown  in  Figure  3.  Thus,  in  order  to  obtain  the 
desired  shape  after  the  nosing  operation,  the  design  of  the  rough-turned  shell 
or  the  preform  prior  to  nosing  is  very  important.  This  problem  is  one  which  costs 
the  shell-forging  industry  thousands  of  dollars,  since  there  is  no  easy  and 
unique  way  to  predict  what  the  original  profile  should  be,  and  several  die  try-outs 
are  necessary  to  establish  the  nosing  process. 

It  has  been  observed  that  certain  sizes  and  shapes  of  shells  elongate 
considerably  during  nosing,  while  in  others  there  is  practically  no  axial  elonga¬ 
tion.  For  those  shells  which  do  not  elongate,  the  problem  is  easier  since  it  can 
be  assumed  that  each  element  of  the  shell  wall  moves  radially  only.  However,  when 
the  shell  elongates  considerably,  the  metal  flow  becomes  complex.  For  shells  with 
uniform  wall  thickness  after  nosing,  Nadai  (Ref.l)  has  outlined  an  approach  for  designing 
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preforms,  based  on  the  strains  in  the  shell  nose.  In  the  following,  Nadai's  approach, 
extended  by  Carlson  (Ref . 4) .  is  used. 


Calculation  of  As-Nosed  Configuration 

Although  as-nosed  configurations  can  be  easily  drawn,  often  numerical 
values  of  the  inside  and  outside  radii  at  a  given  axial  location  is  needed  for 
comparing  measured  values  with  the  theoretical  values.  As  mentioned  earlier, 
the  nosed  portion  of  the  shell  is  described  in  two  different  ways:  (a)  by  two 
radii  and  R^,  as  shown  in  Figure  3(a),  and  (b)  by  a  radius  R^  and  wall 
thickness  h,  as  shown  in  Figure  3(b).  For  the  configuration  given  by  Figure  3(a), 
the  inner  radius  r  and  the  outside  radius  r  at  an  axial  distance  x  is  given  by: 

r.  «  (r  -  h  )  -  R_ (cosa .  -  cosa0)  , 

1  o  o  i  i  2 

r  *  r  -  R,  (cosa  -  cosa,)  ,  ^ 

o  1  o  l 

where  a  ■  sin  1  (B  /R, ) 
o  l  i 

aj^  -  sin  *  (B2/R2) 

-1  Bi  +  x 

a.  -  sin  (  — - -  ) 

1  R! 

-1  B2  +  x 
a  **  sin  (  — - -  ) 

2  R2 

For  the  configuration  described  by  Figure  3(b),  the  outer  radius  is  still  obtained 
as  described  above.  However,  the  inner  surface  is  described  by  (x^,  r^),  the 
coordinates  of  inner  surface  as  below: 


x^  *  (R.^  -  h)  sina^  -  B^ 

r  -  (r  -  h  )  -  (R  -  h) (cosa  -  cosa  ) 
i  o  0  1  o  1 


(2) 
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l. 

r 

!'• 

r 


k  ' 


¥ 


13 

The  above  calculations,  to  determine  the  geometry  of  the  nosed  shell, 
have  been  programmed  on  a  programmable  hand  calculator  (HP-67).  This  program, 
called  NOSE,  can  be  stored  on  a  magnetic  strip  for  easy  reloading  in  the 
calculator.  For  given  geometrical  inputs,  such  as  R^,  B^,  1^,  and  h  or  R^,  B^, 
R^,  B2,  and  l NOSE  calculates  the  coordinates  of  the  inner  and  outer  surfaces 
of  the  nosed  portion  by  incrementing  the  axial  coordinate  by  1  inch.  The  program 
NOSE  and  instructions  for  its  use  are  included  in  Appendix  A. 

In  order  to  illustrate  the  use  of  the  program  NOSE,  calculations  were 
made  for  155  mm  M107  shell.  For  this  purpose,  the  dimensions  were  selected  from 
engineering  drawing  as  given  in  Table  1.  The  calculated  results,  as  displayed  on 
HP-67  calculator,  are  given  in  Table  2. 


Table  1  Dimensions  for  155  mm  M107  Shell 


Outside  radius  at  nose  base,  r 

o 

= 

3 . 10  inch 

Radius  of  ogive,  R^ 

s 

65.50  inch 

Length  of  the  nose  portion,  l ^ 

— 

11.00  inch 

Uniform  wall  thickness,  h 

s 

0.65  inch 

Distance  of  ogive  center  from  nose  base,  B^ 

= 

5.25  inch 

Calculation  of  Preform  Shape 

Figure  4  shows  the  generalized  configuration  of  the  as-nosed  shell,  This 

configuration  will  be  used  in  the  following  design  procedure,  since  it  is 

used  in  defining  large  caliber  shells  where  preform  design  is  crucial.  A  possible 

preform  shape  is  also  shown  in  the  same  figure  with  broken  lines.  The  outside 

surface  of  the  preform  is  straight  and  is  parallel  to  the  axis.  The  inside 

surface  is  made  of  straight-line  segments.  The  length  of  the  preform  is  shorter 

than  the  finish  shape  by  an  amount  e  ,  the  estimated  elongation  of  the  shell 

max 

during  nosing. 
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Table  2  As-Nosed  Profile  of  155  mm  M107  Shell 


Outside  Surface 

Inside 

Surface 

Axial  Distance 

Outer  Radius 

Axial  Distance 

Inner  Radius 

x,  inch 

r ,  inch 

x.,  inch 

r . ,  inch 

X 

i 

0.000 

3.100 

0.000 

2.450 

1.000 

3.012 

0.938 

2.363 

2.000 

2.908 

1.928 

2.260 

3.000 

2.789 

2.918 

2.142 

4.000 

2.654 

3.908 

2.009 

5.000 

2.504 

4.898 

1.860 

6.000 

2.337 

5.888 

1.695 

7.000 

2.155 

6.878 

1.514 

8.000 

1.960 

7.869 

1.318 

9.000 

1.742 

8.859 

1.105 

1C. 000 

1.511 

9.849 

0.877 

11.000 

1.263 

10.839 

0.631 

Fiz  4  Configuration  of  Preform  and  As-Nosed  Shap 
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If  the  shell  does  not  elongate  due  to  nosing,  the  preform  wall  thickness 

h  at  a  distance  x  from  the  nose  base  is  given  as  below: 
o 


rh 


r  cos (a  -  a  ) 
o  o 


* 

•a 


(3) 


If  the  shell  elongates  considerably,  an  element  originally  at  x  goes  to 
x  +  £  after  nosing,  where  £  is  the  axial  displacement  of  the  element.  The  preform 
thickness  in  such  a  case  is  given  by: 


rh 


ho  r  cos (a  -  a  ) 
o  o 


(1  + 


d>f. 


(4) 


The  detailed  derivations  of  the  Equations  (3)  and  (4)  and  a  procedure  for 

estimating  £  from  the  maximum  elongation,  e  ,  are  summarized  in  Appendix  A. 

The  design  of  preform  represents  a  very  important  step  in  nosing  operation 

and  it  will  be  very  useful  to  a  shop  engineer  to  have  the  design  procedure  on  a 

programmable  hand  calculator.  Therefore,  the  above  procedure  for  preform  design 

was  also  programmed  on  an  HP-67.  This  program,  called  PREFORM,  can  be  stored  on  a 

magnetic  strip  for  easy  reloading.  For  given  dimensions  of  the  shell  nose  and 

estimated  elongation  e  ,  beginning  at  the  nose  base,  PREFORM  calculates  preform 

nisx 

wall  thickness  in  increments  of  1  inch.  At  the  end  of  preform  calculation,  if 
the  user  inputs  estimated  values  of  the  coefficient  of  friction  and  the  flow 
stress  of  the  shell  material,  the  program  also  calculates  estimated  nosing  load. 

It  should  be  apparent,  however,  that  the  method  used  for  designing  the  preform 
from  the  known  final  shape  gives  only  approximate  answers.  For  example,  the 
preform  wall  thickness  at  the  tip  cannot  be  determined  using  this  approach,  and 
it  must  be  extrapolated.  The  program  PREFORM  and  instructions  for  its  use  are 
also  included  in  Appendix  A. 
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In  order  to  illustrate  the  use  of  the  program  NOSE,  calculations  were 

made  for  obtaining  the  preform  shape  for  155  mm  M107  shell.  For  this  purpose, 

the  dimensions  given  in  Table  1  were  used,  together  with  an  estimated  value  of 

elongation  e  =  0.500  inch.  The  coefficient  of  friction  at  die-workpiece 
max 

interface  and  the  flow  stress  for  shell  material  were  taken  as  0.10  and  35,000  psi, 
respectively.  Table  3  shows  the  preform  wall  thicknesses  at  various  axial  loca¬ 
tions  from  the  nose  base,  and  the  estimated  nosing  load.  Similar  tables  can  be 
generated  for  other  shells. 


TEMPERATURE  DISTRIBUTION  DUE  TO  PREHEATING  PRIOR  TO  NOSING 


Nosing  of  large  shells  (155  mm  and  above)  is  invariably  done  at  hot- 
working  temperatures.  Hence,  prior  to  nosing,  the  end  portion  of  the  tube  to  be 
formed  is  heated.  Heating  a  portion  of  the  tube  causes  a  temperature  distribution 
along  its  length  prior  to  deformation.  If  the  heating  is  slow  and  gradual,  the 
unheated  portion  of  the  tube  gains  heat  by  conduction.  Consequently,  during  nosing 
the  tube  tends  to  buckle  under  axial  load,  as  shown  in  Figure  2.  Excess  metal  not 
only  accumulates  outside  the  tube,  as  seen  in  Figure  2,  but  also  forms  a  ring  on 
the  inside  near  the  entrance  of  the  nosing  die.  Therefore,  it  is  imperative  that 
heating  is  done  at  a  predetermined  rapid  rate  to  reduce  the  heat  transfer  by 
conduction. 

Because  of  its  numerous  advantages,  induction  heating  is  the  most 
preferred  heating  technique  for  the  nosing  operation.  The  major  advantages  of  the 
induction  heating  are: 

(a)  Uniform  temperature  distribution  all  around  the  shell 

(b)  Precise  control  of  the  final  temperature 

(c)  Rapid  heating  rate,  which  reduces  the  effects  of  scale 

formation  and  heat  conduction.  , 

There  are  many  publications  on  Applied  Mathematics  and  on  Induction 
Heating,  describing  procedures  to  estimate  the  heat  generated  in  induction 
heating  (Ref . 7-14) .  To  avoid  lengthy  derivations,  all  these  procedures  provide 
equations  valid  for  solid  billets  only.  Alternatively,  special-purpose  graphs, 
the  so-called  ’P-Q'  curves,  are  recommended  by  others  for  use  in  the  estimation 


Table  3  Preform  Shape  Prior  to  Nosing  of  155  mm  M107  Shell 


Without  Elongation 

,e  =  0. 

(  max  ) 

With  Elongation 

(emax'  °-5  inch) 

Axial  Distance 
from  Nose  Base, 
inch 

Preform  Wall 
Thickness , 
inch 

Axial  Distance 
from  Nose  Base, 
inch 

Preform  Wall 
Thickness, 
inch 

0.0 

0.650 

0.0 

0.650 

1.0 

0.632 

1.0 

0.632 

2.0 

0.610 

2.0 

0.610 

3.0 

0.585 

3.0 

0.585 

4.0 

0.558 

4.0 

0.558 

5.0 

0.527 

5.0 

0.527 

6.0 

0.492 

6.0 

0.492 

7.0 

0.455 

7.0 

0.455 

8.0 

0.413 

8.0 

0.414 

9.0 

0.369 

9.0 

0.371 

10.0 

0.321 

10.0 

0.354 

11.0 

0.269 

10.5 

0.346* 

Estimated  Nosing  Load, 

P  =  484,355  lbs 

Estimated  Nosing  Load, 

P  *  484,355  lbs 

*  Estimated  by  extrapolation. 
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of  heat  generated  in  tubular  charges (Ref • 13-17) .  However,  the  existing  P-Q  curves 
are  valid  only  for  non-magnetic  materials  and  additional  curves  need  to  be 
generated  for  magnetic  materials. 

A  brief  look  at  the  references  cited  above  will  indicate  that  most  of 
the  analyses  on  induction  heating  have  been  conducted  in  a  period  when  use  of 
computers  for  scientific  analyses  was  fairly  unknown.  Therefore,  simplified 
procedures  by  approximate  methods  were  necessary  to  avoid  time-consuming 
calculations.  Furthermore,  most  of  these  analyses  were  aimed  at  the  design  of 
induction  coils  for  specific  applications (Ref . 15-17) .  Hence,  the  given  equations 
are  good  only  to  estimate  the  total  heat  generated  in  a  given  time  interval. 

In  recent  years,  some  companies  specializing  in  induction  heating  have 
improved  existing  analytical  procedures  and  developed  computer  programs  for 
advanced  analysis  (Ref . 18-20) .  By  use  of  these  computer  programs,  temperature  distri¬ 
butions  in  slab,  cylindrical,  and  tubular  charges  made  up  of  magnetic  and 
non-magnetic  materials  can  be  estimated  under  a  variety  of  heating  conditions. 
However,  neither  the  analytical  procedure  nor  the  computer  programs  are  readily 
available  for  public  use.  In  addition,  heating  of  tubes  prior  to  nosing  poses 
another  problem.  In  this  operation,  only  a  part  of  the  tube  is  heated.  During 
the  heating  and  the  subsequent  transfer  of  the  tube  to  the  nosing  press,  heat  is 
conducted  to  the  unheated  portion  of  the  tube.  Thus,  the  heat  transfer  causes  a 
temperature  gradient  along  the  axis  of  the  tube.  The  nature  of  this  temperature 
gradient  is  very  critical  to  avoid  the  buckling  of  the  tube  during  nosing.  The 
existing  general  equations  (Ref . 21-22)  are,  therefore,  not  applicable  to  partially 
heated  tubes.  Hence,  in  order  to  estimate  the  temperature  distribution  in  a 
partially  heated  tube,  special-purpose  equations  were  derived.  The  derivation 
procedure  and  the  derived  equations  are  given  in  Appendix  B. 

Heat  Transfer  During  Induction  Heating 

During  induction  heating  of  a  tube  end,  a  major  portion  of  the  heat 
generated  remains  in  the  charge,  i.e.,  the  tube;  a  part  of  the  heat  flows  into 
the  unheated  cold  section  of  the  tube  by  conduction;  an  additional  fraction  of 
the  heat  generated  is  lost  to  the  environment  by  convection  and  radiation.  In 
order  to  predict  the  local  temperatures,  it  is  necessary  to  consider  the  heat 
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equalization  due  to  heat  transfer.  Because  of  the  complexity  of  the  phenomenon, 
it  is  appropriate  to  use  the  finite-difference  method  for  solving  the  heat 
transfer  equations.  During  a  small  time  interval,  At,  it  is  assumed  that  the 
heat  generation  takes  place  instantaneously  at  the  beginning  of  the  interval  At. 
Using  these  calculated  values  of  temperature  as  temperatures  before  At,  the  heat 
flow  is  analyzed  and  the  temperature  field,  which  exists  after  the  heat  loss  by 
conduction,  convection,  and  radiation  during  the  same  interval  At,  is  determined. 
The  repetition  of  these  two  steps  simulates  numerically  the  heat  generation  and 
transfer  simultaneously  and  gives  the  temperature  distribution  as  a  function  of 
time.  The  difference  equations  used  in  the  finite-difference  heat  transfer 
analysis,  under  different  boundary  conditions  of  the  shell-nosing  operation,  are 
also  given  in  Appendix  B. 

Heat  Transfer  after  Induction  Heating 

After  the  heating  of  the  tube,  temperature  equalization  takes  place 
during  the  interval  between  heating  and  the  actual  forming  operation.  During  this 
period,  no  heat  is  generated,  but  heait  is  conducted  to  the  cooler  portion  of  the 
tubes  and  is  convected  and  radiated  to  the  atmosphere.  Temperature  fields  which 
exist  during  and  after  this  cooling  period  are  generated  using  the  same  procedure 
and  the  equations  used  before,  but  excluding  the  heat  generation  aspect. 

Prediction  of  Temperature  Fields  in 
Inductively  Heated  Tubular  Components 

In  order  to  effectively  apply  these  principles  of  heat  generation  and 
transfer  to  practical  problems,  the  entire  procedure  is  computerized.  The  set 
of  computer  programs  performing  this  analysis  is  called  INHEAT.  Figure  5  shows 
the  functional  flow  chart  of  INHEAT.  As  seen  in  Figure  5,  a  complete  simulation 
of  the  induction  heating  and  subsequent  temperature  equalization  is  undertaken 
to  determine  the  non-steady  state  temperature  distributions  in  tubes  at  various 
time  intervals.  This  process  analysis  is  applicable  to  both  magnetic  and  non¬ 
magnetic  materials. 
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Fig  5  Functional  Flow  Chart  of  the  Computer  Program  INHEAT 
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For  the  purpose  of  heat  transfer  analysis  by  finite-difference  method, 
a  portion  of  the  vertical  cross  section  of  the  tube  is  considered  as  to  be 
divided  into  a  number  of  trapezoidal  grid  elements,  as  seen  in  Figure  6.  The 
maximum  value  of  At  for  the  proposed  grid  pattern  is  determined  using  the 
stability  criterion,  given  in  Appendix  B.  For  simulation,  it  is  assumed  that  the 
temperature  rise  due  to  induction  heating  takes  place  instantaneously  during  the 
time  interval  At,  followed  by  heat  transfer  during  the  same  time  interval. 

At  the  beginning,  the  temperature  in  each  grid  point  is  set  to  the  same 
value  as  the  preheat  temperature  of  the  tube,  or  the  room  temperature  if  the  tube 
is  not  preheated.  Temperature  rise  due  to  induction  heating  during  a  time  interval 
At  is  calculated  from  the  process  variables  and  added  to  each  grid  within  the 
heating  zone  of  the  tube.  Using  these  calculated  values  of  temperature,  the  heat 
flow  is  analyzed  and  the  temperature  field  which  exists  after  heat  transfer  in  the 
time  interval  At  is  determined.  By  repeating  the  sequence  of  heat  generation  and 
heat  transfer,  the  non-steady  state  temperature  fields  in  the  tube  at  various  time 
intervals  are  determined. 

During  this  analysis,  all  the  temperature  dependent  material  properties, 
such  as  heat  capacity,  thermal  conductivity,  electrical  resistivity  and  magnetic 
permeability  are  estimated  separately  for  each  grid  in  the  system  using  its 
instantaneous  temperature. 

Following  induction  heating,  temperature  equalization  takes  place  during 
the  transfer  of  the  billet  to  the  press.  Thus,  the  temperature  distribution  prior 
to  nosing  is  the  one  after  temperature  equalization.  Hence,  to  obtain  the  tempera¬ 
ture  distribution  prior  to  nosing,  the  same  heat  transfer  analysis,  excluding  the 
heat  generation  aspect,  is  performed  for  the  duration  of  billet  transfer.  After 
the  simulation,  temperature  distributions  at  the  end  of  the  induction  heating  and 
prior  to  nosing  are  printed  in  a  tabular  form  as  shown  in  Figure  B-2  of  Appendix  B, 
attached  to  this  report.  Intermediate  results,  if  desired,  may  also  be  obtained 
by  assigning  appropriate  values  to  program  control  parameters,  as  described  in 
Appendix  B. 
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Experimental  Evaluation 

■* 

To  assess  the  applicability  and  the  accuracy  of  the  computer  programs 
INHEAT,  the  following  special-purpose  experiments  were  conducted.  As  illustrated 
in  Figure  7,  four  thermocouples  were  welded  on  the  outer  surfaces  of  each  of  two 
2-inch  OD,  1.5-inch  ID  and  8-inch  long  tubes.  One  of  the  tubes  was  made  of 
stainless  steel  type  304.  This  tube  was  totally  annealed  (15  minutes  at  1050  C 
and  rapidly  cooled  in  air  blast)  prior  to  induction  heating  to  eliminate  the 
effect  of  cold  work  on  its  magnetic  permeability.  The  other  tube  was  made  of 
AISI  1045  medium  carbon  steel.  An  induction  coil,  with  14  turns  over  a  length 
of  four  inches,  was  specially  made  using  0.25-inch  OD  copper  tube.  For 
insulation,  the  entire  length  of  the  tube  was  wrapped  with  fiberglass  sheets. 

To  eliminate  the  excessive  heating  of  the  induction  coil  during  induction 
heating,  water  was  continuously  circulated  through  the  tube  using  a  separate 
pump. 

Figure  8  illustrates  the  experimental  set  up  with  one  of  the  steel 
tubes  in  position  in  the  induction  coil.  The  terminals  of  the  thermocouples 
are  connected  to  two,  two-channel  Honeywell  Electronic  19  recorders.  The  • 

power  to  the  induction  coil  was  supplied  by  an  Ajax  Magnathermic  Corporation's 
induction  machine.  Model  23HT.  This  machine  is  equipped  with  a  motor  generator 
set.  The  output  voltage  of  the  generator  can  be  adjusted  between  0  and  800 
volts;  consequently,  the  output  current  and  the  output  power  will  vary  between 
0  and  125  amps  and  0  and  100  KVA,  respectively.  To  improve  the  power  factor, 
the  generator  output  circuit  is  equipped  with  an  adjustable  capacitor  set. 

The  output  circuit  is  also  equipped  with  an  adjustable  transformer  to  match  the 
load  voltage  of  the  coil  to  the  generator  output  voltage.  The  output  frequency 
of  the  generator  is  fixed  at  10  KHz.  The  equivalent  heating  circuit  used  in 
this  equipment  is  schematically  represented  in  Figure  9. 

Prior  to  the  actual  experiment,  several  preliminary  tests  were  run 
(a)  to  improve  the  power  factor  of  the  generator  output  by  adjusting  the 
capacitor  in  the  circuit,  and  (b)  to  match  the  load  voltage  to  the  output 
voltage  of  the  generator  by  varying  the  transformer  ratio.  During  the 

*  . 

individual  heating  of  both  steel  tubes,  output  from  all  the  four  thermocouples 

were  recorded  on  Honeywell  recorders.  The  current  in  the  induction  coil  circuit  j 

i 
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Experimental  Set  Up  for  Induction  Heating  of  Tubes 


Equivalent 
load  resistance 
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could  not  be  measured  due  to  the  non-availability  of  suitable  meter  to  measure 
1000  amps  at  10  KHz.  However,  the  voltage  across  the  coil  was  measured  using 
a  Hewlett-Packard  RMS  voltmeter  capable  of  functioning  between  50  Hz  and  3  MHz. 
The  output  voltage  of  the  generator  was  adjusted  such  that  the  tubes  will  be 
heated  from  room  temperature  to  about  1000  C  in  about  20  seconds. 


Computer-Aided  Estimation  of  Temperatures  for  Experimental  Conditions 


Since  the  current  in  the  induction  coil  circuit  could  not  be  measured 
directly,  it  was  estimated  by  the  following  procedure: 

From  the  skin  depth  of  the  tube,  the  load  resistance  was  calculated 
using  the  formula: 


p  ff  D 
5  J l 


where 


and 


o  is  the  electrical  resistivity  of  the  material 
D  is  the  outer  diameter  of  the  tube 


i3  the  length  of  the  induction  coil 
6  is  the  skin  depth  and  is  given  by 


_1 
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Using  the  transformer  ratios,  the  equivalent  load  resistance  in  the 

primary.  Figure  9,  was  estimated  from  R,  given  by  Equation  (5).  The  inductance 

X£  was  calculated  from  the  capacitor  in  the  primary  circuit.  Since  the  circuit 

was  tuned  up,  it  was  assumed  that  the  inductance  of  the  adjustable  transformer 

primary,  Xp  in  Figure  9,  is  equal  to  Xc  and  its  resistance  is  negligible. 

From  the  known  primary  input  voltage  E  and  the  resistance  R^  and  inductance  X  , 

P  1>  P 

the  current  in  the  primary  was  obtained  applying  Ohm's  law.  The  current  in 

the  secondary,  I  ,  is  the  product  of  the  primary  current  I_,  and  the  transformer 
c  l 

ratio. 

Also,  since  the  secoftdary  voltage  was  measured  separately,  the  secondary 
current  I  was  estimated  independently  following  the  procedure  described  in 
References  15  and  16.  Both  the  estimated  currents  were  close  within  acceptable 
engineering  accuracy. 
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Using  this  estimated  current  value  in  INHEAT,  temperatures  at  various 
points  in  the  tube  were  estimated.  The  formula  given  in  Appendix  B  to  estimate 
the  temperature  rise  due  to  induction  heating  is  for  ideal  heating  condition. 

Because  the  practical  conditions  differ  considerably  from  the  ideal  situation, 
an  efficiency  factor  was  introduced  in  Equation  (B-lla)  of  Appendix  B  in 
estimating  the  temperatures.  This  efficiency  factor,  which  depends  upon  (a) 
the  coil  design,  (b)  cross-sectional  geometry  of  the  copper  tube  forming  the 
coil,  (c)  air  gap  between  the  coil  and  the  tube  heated,  and  (d)  physical 
properties  of  the  material  being  heated,  was  estimated  as  50  percent  in  the 
present  tests.  In  a  practical  condition,  the  efficiency  factor  can  be  estimated 
by  running  a  preliminary  trial,  measuring  the  temperature  at  a  specific  point 
on  the  tube  and  comparing  the  results  with  the  estimated  temperature. 

Comparison  of  Predictions  with  Experimental  Results 

For  comparison  with  the  actually  measured  temperatures,  the  estimated 
temperatures  at  different  locations  of  the  thermocouples  were  plotted  along  with 
the  measured  temperatures.  Figure  10  shows  the  temperature  distributions  in  type 
304  stainless  steel.  Estimated  temperatures  in  thermocouples  1,  2  and  3  (refer 
to  Figure  7  for  the  locations  of  the  thermocouples)  are  very  close  to  each  other 
at  all  times.  Measured  temperatures  in  thermocouples  1  and  2  are  similarly 
distributed  as  the  estimated  temperatures,  except  for  a  slow  response  during  the 
early  stages  of  heating.  This  slow  response  is  attributed  to  the  inertia  effect 
of  the  mechanical  type  recorder  used  in  this  experiment.  With  the  use  of  optical 
or  some  other  quick  response  recorders,  this  discrepancy  could  be  eliminated. 

The  measured  temperature  in  thermocouple  3  is  very  much  lower  than  the 
estimated  temperature.  In  the  estimation  procedure,  temperature  rise  due  to 
induction  heating  is  assumed  to  be  uniform  over  the  entire  length  of  the  induction 
coil.  However,  when  the  length  of  the  induction  coil  is  shorter  than  the  workpiece, 
heat  generated  closer  to  the  end  of  the  induction  coil  is  less  than  predicted,  due 
to  end  effects.  Due  to  this  well-known  phenomenon,  a  provision  has  already  been 
incorporated  in  the  computer  programs  to  account  for  the  variations  in  heat 
generation  along  the  axis  of  the  coil;  but  due  to  the  lack  of  any  reliable  data, 
this  particular  feature  could  not  be  used  in  the  present  analysis.  As  a  result. 
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the  estimated  temperature  distribution  at  the  location  of  thermocouple  3  is  much 
higher  than  the  measured  distribution.  In  order  to  determine  the  true  distribu¬ 
tion  of  heat  generation  along  the  axis  of  the  coil,  it  is  necessary  to  run  more 
trials  with  a  number  of  thermocouples  along  the  entire  length  being  heated. 

A  similar  trend  of  deviations  exhibited  by  thermocouple  4  (refer  to  Figure  10) 
is  purely  a  consequence  of  the  difference  between  the  measured  and  estimated 
temperature  distributions  from  thermocouple  3. 

Figure  11  shows  the  temperature  distributions  in  AISI  1045  steel. 

Similar  to  the  results  in  Figure  10,  the  retarded  response  of  the  measured 
temperatures  may  be  attributed  to  the  inertia  effect  of  the  mechanical  recording 
device.  Unlike  in  the  stainless  steel,  the  measured  temperature  distribution 
from  thermocouple  1  is  very  much  lower  than  the  estimated  distribution.  Although 
the  above  discrepancy  may  be  attributed  to  various  factors  associated  with  the 
assumptions  of  the  mathematical  model  and  with  the  conditions  of  the  experiments. 
However,  the  main  factor  responsible  for  poor  correlation  between  the  predicted 
and  the  experimental  results  near  the  lower  end  of  the  coil,  Figure  12,  is  the 
well-known  "end  effect".  In  the  mathematical  modeling  of  the  problem,  this  end 
effect  was  assumed  to  be  negligible  for  the  sake  of  simplicity  and  a  one-dimen¬ 
sional  analysis  was  developed.  A  rigorous  two-dimensional  analytical  approach, 
considering  the  end  effect,  may  be  very  complicated  and  may  or  may  not  signifi¬ 
cantly  improve  the  end  results.  Hence,  it  is  believed  that  a  combined  analytical- 
empirical  approach  may  be  easier  and  more  suitable  for  the  present  application. 

Because  of  the  end  effects,  the  amount  of  heat  generated  per  unit  volume 
near  the  lower  end  of  the  coil  is  smaller  than  that  around  the  mid  length  of  the 
coil.  Thus,  the  generation  of  induction  heat,  due  to  a  uniformly  spaced  coil,  is 
not  uniform  along  the  length  of  the  coil  and  will  be  as  shown  in  Figure  12.  This 
pattern  of  heat  generation  may  depend  upon  several  factors,  such  as  (a)  size  of 
the  tube,  (b)  magnetic  properties  of  the  material  heated,  (c)  the  frequency  of  the 
heating  current,  (d)  heating  rate,  and  more  importantly,  (e)  the  coil  design. 
Neither  analytical  nor  empirical  results  considering  the  end  effects  are  available 
in  any  reviewed  technical  publications.  Hence,  it  is  necessary  to  conduct  experi¬ 
ments  with  105  mm,  155  mm,  175  mm  and  8-inch  diameter  shell  preforms  to  generate 
the  heat  generation  profiles.  At  this  time,  it  is  expected  that  in  a  given 
installation  (i.e.,  for  a  particular  coil  design),  the  size  of  the  tube  on  the 
heat  generation  pattern  will  have  the  maximum  effect  on  heat  generation  and  the 


Fig  10  Measured  and  Computed  Temperature  Distributions  in  Stainless  Steel  Tube 


Fig  11  Measured  and  Computed  Temperature  Distributions  in  1045  Steel  Tube 


other  process  variables,  listed  above,  will  have  negligible  effects.  By  conducting 
a  planned  set  of  experiments,  the  least  important  variables,  influencing  heat 
generation  in  the  shell  wall,  can  be  identified  and  eliminated  from  the  list. 

Also,  the  heat  generation  patterns  under  various  conditions  of  the  effective 
variables  can  be  established.  Once  these  profiles  are  determined,  the  computer 
program  INDHET  can  use  these  heat  generation  patterns  as  inputs  and  predict  the 
results  more  accurately  than  the  predictions  now. 

ANALYSIS  AND  PREDICTION  OF  METAL  FLOW  IN  NOSING 

In  nosing  of  shells,  the  shape  of  the  outer  surface  of  the  shell  follows 
the  die  profile  while  the  inner  surface  is  not  supported,  as  shown  in  Figure  13. 
Thus,  during  the  nosing  operation,  the  part  of  the  tube  or  shell  within  the  die 
zone  can  move  either  radially  inward,  or  in  the  direction  of  the  tangent  to  the 
die  profile,  or  in  both  directions,  simultaneously.  Under  plastic  state,  metals 
flow  in  the  direction  of  least  resistance;  therefore,  the  metal  flow  in  cold  nosing 
is  affected  by  the  frictional  restraint  at  the  die-workpiece  interface  and  by  the 
work-hardening  characteristics  of  the  material  being  deformed.  In  hot  nosing,  the 
metal  flow  depends  upon  the  friction,  the  temperature  of  the  workpiece  and  the 
speed  of  the  die.  In  general,  wall  thickening  increases  with  increasing  friction, 
and  in  fact,  under  severe  friction  conditions,  the  nosed  tube  or  shell  may  shrink 
in  length  instead  of  elongating. 

Under  this  program,  a  mathematical  model  for  predicting  metal  flow  in 
nosing  of  shells  was  developed.  This  model  considers  preform  with  uniform  wall 
thickness  and  utilizes  Hill's  general  method  of  analysis  for  metalworking 
processes  (Ref . 23-24) .  Thus,  the  model  is  exactly  valid  for  nosing  of  shells  up  to 
105-mm  size,  which  are  usually  nosed  from  preforms  with  uniform  walls.  For 
larger  shells,  where  the  preform  wall  thickness  is  not  uniform,  this  analysis  can 
be  applied  only  approximately. 


Fig  13  Configuration  of  Shell  During  Nosing 
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Analysis 


The  first  step  in  using  Hill's  method  consists  of  choosing  a  class  of 

kinematically-admissible  velocity  fields  (with  h(a),  the  unknown  thickness 

function  in  a  (see  Figure  13)  from  which  the  best  approximation  will  eventually 

be  selected.  The  associated  stress  distribution  is  determined  using  the  flow 

rule  and  the  yield  criterion,  which  will  generally  not  satisfy  all  the  statical 

requirements.  Therefore,  a  second  kinematically-admissible  velocity  field,  called 

the  orthogonalizing  field,  and  which  can  be  similar  to  the  first  velocity  field, 

is  selected  and  the  converse  to  the  virtual  work-rate  principle  is  applied.  This 

results  into  a  system  of  equilibrium  equations  and  boundary  conditions,  suited  to 

the  particular  approximating  velocity  field  and  uniquely  determining  its  best 

member.  The  details  of  Hill's  method  and  the  analysis  of  metal  flow  in  nosing  of 

shells  using  this  method  are  included  in  Appendix  C. 

As  given  by  Equation  (C-15)  in  Appendix  C,  the  thickness  function,  h(a), 

of  shell  wall  is  defined  by  a  second  order  ordinary  differential  equation.  Using 

the  boundary  conditions  given  by  Equation  (C-16)  in  Appendix  C,  Equation  (C-15)  can 

be  easily  solved  using  a  numerical  technique.  However,  since  the  boundary  conditions 

are  available  at  each  end  of  the  interval,  some  type  of  interative  method  is  required 

for  numerical  integration  of  Equation  (C-15) .  For  this  purpose,  initially  two 

guesses  of  h'  (aQ)  are  made  and  the  Equation  (C-15)  is  solved  by  a  fifth  order 

Runge-Kutta  Method.  Normally,  these  solutions  will  not  satisfy  the  second  of  the 

boundary  conditions,  Equation  (C-16).  Therefore,  a  third  guess  value  of  h'  (a  ) 

o 

is  determined  using  the  first  two  solutions  by  linear  extrapolation  such  that  the 

projected  error  in  h'(a^)  is  zero,  and  the  entire  integration  procedure  is 

repeated.  This  last  step  is  repeated  with  last  two  solutions  until  the  error  in 

h'  (a^)  is  within  specified  error  bounds.  Further,  the  Equation  (C-15)  is  singular 

when  a  =0,  since  r'  (a  =  a  *  0)  *  0.  In  this  particular  case,  a  solution  can 
o  o 

be  obtained  by  taking  uq  as  a  small  positive  quantity  instead  of  zero.  The  error 

caused  by  this  approximation  is  relatively  small. 

Based  on  the  above  analysis,  a  computer  program,  named  NOSFLW,  was 

developed  to  predict  metal  flow  during  nosing  of  shells.  The  program  NOSFLW  is 

coded  in  FORTRAN  IV  and  requires  approximatelv  25,000o  words  of  memory  space  in  a 

o 
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CDC  CYBER  70  computer.  A  typical  run  requires  approximately  7.0  system  seconds, 
including  5.5  seconds  on  the  central  processor.  All  the  input  data  to  the 
computer  program  NOSFLW  are  transferred  through  READ  statements.  These  include 
variables  defining  the  geometries  of  the  preform  and  nosed  shell,  the  friction 
shear  factor  at  the  die-workpiece  interface,  the  velocity  of  the  die,  the  flow 
stress  of  the  workpiece  material  and  two  guess  values  of  the  slope  of  the  inner 
profile  at  the  base  of  the  nose.  As  output,  the  computer  program  prints  the 
coordinates  of  the  inner  and  outer  profiles  of  the  nose,  the  elongation  due  to 
nosing,  the  maximum  thickening  of  the  wall  and  the  estimated  nosing  load.  The 
computer  printout  also  gives  the  velocity  field  in  the  deformation  zone,  i.e., 
the  velocities,  the  strains  and  strain  rates  as  functions  of  locations  along  the 
nosed  portion.  A  functional  flow  chart  of  the  computer  program  NOSFLW  is  given 
in  Figure  14. 


Parametric  Study 


In  order  to  illustrate  the  application  of  the  analysis  and  the 
associated  computer  program,  predictions  were  made  for  cold  nosing  of  105-mm 
Ml  shell  from  a  tubular  preform  with  uniform  wall  thickness.  For  this  purpose, 
the  values  of  the  input  variables  to  the  computer  program  NOSFLW  were  selected 
from  the  engineering  drawings  (after  including  the  finish  machining  allowances) 
as  given  in  Table  4. 


Table  4  Dimensions  for  105-ram  Ml  Shell 


Outside  diameter  at  the  nose  base 
Radius  of  curvature  of  the  nose  ogive 
Initial  wall  thickness  of  the  preform 
Axial  length  of  the  nosed  portion 
Distance  of  the  ogive  center  from  nose  base 


=  108  mm 
=  650  mm 
*  12  mm 
=  150  mm 
=  25  mm 


39 


The  friction  shear  factor  at  the  die-workpiece  interface  was  estimated 
to  be  0.05.  The  computer  program  predicts  a  maximum  wall  thickening  of  68.34 
percent  and  an  elongation  of  1.7  percent  after  nosing  is  completed.  The  actual 
preform  geometry  and  the  predicted  shell  geometry  after  nosing  is  shown  in 
Figure  15.  Although  no  experimental  data  were  available  for  comparison,  the 
nosed  shape  in  Figure  15  is  in  reasonably  good  agreement  with  the  desired  shape 
after  nosing  for  105-mm  Ml  shell. 

As  stated  earlier,  depending  upon  the  frictional  restraint  at  the 
tool-workpie^e  interface,  the  shell  length  may  increase  or  decrease  as  it 
penetrates  in  the  nosing  die.  As  seen  in  Figure  16,  at  low  values  of  the 
friction  shear  factor,  the  shell  elongates  due  to  nosing,  whereas,  at  high 
values  of  friction,  the  shell  elongates  in  the  initial  stages  of  nosing  and  then 
starts  to  shrink  in  length  as  penetration  into  the  die  is  increased.  This 
predicted  trend  is  in  agreement  with  experimental  observations  reported  in  the 
literature  (Ref . 6) . 

During  nosing,  the  shell-wall  thickness  invariably  increases  with 
increasing  penetration  in  the  die  and  with  increasing  frictional  restraint  at  the 
die-workpiece  interface.  Figure  17  illustrates  this  trend  as  predicted  by  the 
present  analysis  and  the  associated  computer  programs. 

As  a  by-product,  the  present  analysis  also  predicts  the  load  required 
during  nosing  of  shells.  Figure  18  shows  predicted  values  of  the  nosing  load 
normalized  with  respect  to  area  of  cross  section  and  the  flow  stress  at  the  nose 
base  for  various  values  of  the  friction  shear  factor  selected  in  the  present 
study.  The  trend  is  in  agreement  with  experimental  observations  (Ref. 4).  This 
analysis  is  also  capable  of  predicting  the  limiting  conditions  of  the  process. 

As  seen  in  Figure  18,  a  complete  nosing  operation  will  not  be  possible  with  a 
friction  shear  factor  larger  than  0.07  due  to  local  yielding  at  the  nose  base. 

For  this  reason,  good  lubrication  is  essential  in  nosing  of  shells,  apart  from 
obtaining  a  desired  geometry  after  nosing. 


- Preform 


As  nosed 
lm«003) 


Ft*  15  Preform  Shape  and  A.s-'Sosed  Shape  Predicted  by 

the  Computer  Program  NOSFLW  for  105  mm  Ml  Shell 

(All  Dimensions  in  mm) 
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Penetration  in  the  Die,  mm 


Fig  16  Elongation  in  Shell  Length  Due  to  Nosing 
at  Various  Friction  Values 
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Increase  in  Wall  Thickness  Percent 
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Penetration  in  the  Die,  mm 


Fig  17  Increase  in  Wall  Thickness  Due  to 
Nosing  at  Various  Friction  Values 


Penafrotion  in  the  Die, mm 


Fig  18  Normalized  Load  During  Nosing  at  Various  Values 

of  Friction  Shear  Factor  (p  =  Nosing  Load,  A  and 
a  are  the  area  of  cross  section  and  the  flow 
stress  at  the  nose  base) 
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ANALYSIS  OF  STRESSES  AND  MATERIAL  FAILURE 
IN  NOSING  OF  SHELLS 

Nadai(Ref.l)  analyzed  the  stresses  in  nosing  of  shells.  The  analysis  assumed 
that  the  tube  material  is  perfectly  plastic  and  that  the  wall  thickness  does  not 
change  appreciably  during  the  nosing  process.  Onat  and  Prager  (Ref . 2)  extended  Nadai's 
work  and  included  the  changes  in  the  shell-wall  thickness.  However,  both 
analyses  assume  a  preform  of  uniform  wall  thickness  and  calculate  the  state  of 
stresses  at  the  end  of  the  nosing  operation.  Therefore,  in  order  to  generate  a 
load-stroke  curve  and  treat  preforms  with  nonuniform  wail  thickness,  a  step-by- 
step  approach  of  stress  analysis  is  described  in  the  following.  In  addition,  the 
following  analysis  is  valid  for  real  materials  since  it  considers  the  flow  stress 
of  the  deforming  material  as  a  function  of  strain,  strain  rate  and  temperature. 

A  schematic  diagram  of  nosing  of  shells  is  shown  in  Figure  19.  The 
length  of  the  preform  is  divided  into  a  number  of  segments.  For  nonuniform-walled 
preforms,  if  the  number  of  segments  is  sufficiently  large,  the  variation  in 
thickness  along  the  length  of  the  segment  can  be  ignored,  and  the  average  thickness 
of  the  segment  can  be  taken  as  its  uniform  thickness.  Using  this  simplification, 

Onat  and  Prager 's  analysis  can  be  applied  to  each  element.  If  the  preform  has 
uniform  wall  thickness,  this  simplication  is  not  necessary.  However,  for  the 
sake  of  generality,  the  analysis  was  developed  for  nosing  with  nonuniform-walled 
preforms . 

The  nosing  operation  can  be  simulated  by  moving  the  die  over  one  segment 
of  the  preform  at  a  time.  Thus,  when  a  segment  is  inside  the  die,  it  is  deformed 
plastically  and  takes  a  new  average  thickness  given  by  (3),  if  the  shell  does  not 
elongate,  or  by  Equation  (4),  if  elongation  due  to  nosing  is  considered.  During 
the  next  step,  this  segment  is  pushed  further  inside  the  die  and  a  new  element  of 
the  preform  is  deformed  for  the  first  time.  At  each  step,  the  stresses  and  the 
thickness  of  each  element  are  calculated.  Thus,  the  nosing  process  is  simulated 
in  a  finite  number  of  discrete  steps. 
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Analysis  of  Stresses 

Let  us  consider  the  analysis  of  stresses  for  element  (1),  Figure  19. 

The  analysis  of  stresses  for  the  subsequent  elements  will  be  the  same,  except  for 

the  boundary  conditions  which  are  to  be  determined  from  the  known  stress  distribution 

of  the  preceding  element.  The  surface  traction  exerted  by  the  die  on  the  nose  is  the 

normal  pressure,  p,  and  the  friction  shear  stress,  f  =  yp,  where  y  is  the 

coefficient  of  friction.  The  meridianal  stress  and  the  circumferential  stress 

produced  by  these  surface  tractions  in  the  thin-walled  nose  far  exceed  p  and  f 

in  absolute  value.  Therefore,  the  state  of  stress  in  the  nose  wall  can  be 

treated  as  plane  with  the  circumferential  stress  o  and  the  meridianal  stress 

c 

o  as  the  principal  stresses.  Since  these  principal  stresses  are  compressive, 
m 

and  since  a  can  be  expected  to  exceed  a  in  absolute  value,  according  to 
cm' 

Tresca's  yield  criterion, 

oc  -  -  o  (6) 

where  a  is  the  material  flow  stress,  and  it  is  a  function  of  the  strain,  strain 
rate  and  temperature.  Equilibrium  of  forces  in  the  meridianal  and  the  circum¬ 
ferential  directions,  together  with  Equation  (9),  yield  the  following: 


do 
_ m 

da 


+  y  a  =  a 


sina  9-  y  cosa 
cosa  -  a/R 


(7) 


Integration  above  equation  gives  the  distribution  of  the  meridianal  stress  a  as 

m 

below: 
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da 

(8) 


where  o'  is  the  meridianal  stress  at  a  =  a^,  the  tip  of  the  element  under  considera¬ 
tion.  For  element  (1),  o  =  0  at  a  =  a, .  With  this  boundary  condition,  the 

mi 

meridianal  stress  a  can  be  found  as  a  function  of  a  for  element  (1).  For  the 
m 

subsequent  elements,  the  boundary  condition  at  the  leading  end  of  the  element 
is  given  from  the  known  value  of  from  the  analysis  of  stresses  in  the  preceding 


element.  Thus,  because  the  value  of  a  at  the  boundary  of  elements  (1)  and  (2) 

m 

is  known  from  the  analysis  of  stresses  in  element  (1),  Equation  (8)  can  be 
solved  to  determine  stresses  in  element  (2).  Similar  procedure  is  applied  to 
subsequent  elements.  A  detailed  derivation  of  this  stress  analysis  is  given  in 
Appendix  D. 

In  the  calculations,  the  flow  stress  of  each  element  considered  as  a 
function  of  the  strain,  strain  rate,  and  temperature.  The  strain  in  an  element 
is  the  cumulative  strain.  The  strain  rate  is  calculated  from  a  simple  velocity 
field  given  in  Appendix  C.  The  temperature  of  an  element  will  depend  upon  the 
heat  generated  due  to  plastic  deformation  and  friction,  and  the  heat  conduction 
to  the  colder  dies.  The  expressions  for  the  strain,  strain  rate  and  the  tempera¬ 
ture  of  an  element  inside  the  die  are  also  included  in  Appendix  D. 

Material  Failure 

In  practical  nosing  operations,  material  failure  sometimes  occur  and  it 

has  been  reported  from  various  shell-forging  plants  (Ref. 6).  In  hot  nosing,  the 

critical  load  is  a  function  of  the  temperature  distribution  along  the  nose  of  the 

shell  as  well  as  the  wall  thickness  to  diameter  ratio.  Thus,  if  the  temperature 

distribution  is  not  correct,  local  bulging  at  the  base  of  the  nose  can  occur  due 

to  improper  metal  flow.  This  local  swelling  of  the  wall  thickness  can  also  occur 

in  the  curved  portion  of  the  nose.  To  avoid  this  type  of  material  failure,  the 

condition  to  be  satisfied  is  that  the  meridian  stress  a  at  the  base  of  the  nose 

m 

should  not  reach  the  yield  stress  in  pure  compression  at  the  nose  base.  Thus,  to 
avoid  local  bulging  in  hot  nosing, 

a  <  a  at  a  =  a  .  (9) 

m  o 

In  cold  nosing,  unless  the  shell  is  properly  designed,  it  may  buckle 
under  the  forces  necessary  for  nosing.  In  most  cases,  the  shell  wall  is  so  thick 
that  the  proportional  limit  of  the  material  will  be  reached  before  the  equilibrium 
of  the  shell  becomes  unstable,  and  the  local  buckling  will  occur.  The  elastic 
buckling  of  cylindrical  shells  is  treated  by  Timoshenko.  The  critical  stress 
necessary  to  cause  elastic  buckling  is  given  by  the  following  formula: 
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cr 


E  h _ 

r  /3(1  -  v  2) 
o 


(13) 


A 


where  E  =  modulus  of  elasticity 

h  =  thickness  of  shell 

r  =  outside  radius  of  shell 
o 

u  =  poisson's  ratio. 

This  formula  is  valid  when  the  critical  stress  is  below  the  proportional  limit  of 
the  material.  Thus,  at  each  step  of  simulation,  both  local  bulging  and  elastic 
buckling  is  investigated. 

Computer  Simulation 


In  order  to  facilitate  the  application  of  the  stress  and  the  material 
failure  analyses,  a  computer  program,  called  NOSTRS,  was  developed  to  simulate 
the  nosing  process.  This  computer  program  is  capable  of  treating  both  cold  and 
hot  nosing  of  shells.  A  function  flow  chart  of  the  computer  program  NOSTRS  is 
given  in  Figure  20.  The  computer  program  NOSTRS  is  coded  in  FORTRAN  IV  and 
requires  approximately  35,000g  words  of  memory  space  in  a  CDC  CYBER-70  computer. 

A  typical  run  requires  approximately  16  seconds,  including  9.7  seconds  on  the 
central  processor.  All  the  input  data  to  the  computer  program  NOSTRS  are  trans¬ 
ferred  through  READ  statements.  These  include  variables  defining  the  geometry  of 
the  nosed  shell,  the  friction  coefficient  at  the  die-workpiece  interface,  the  die 
velocity,  the  type  of  preform  (uniform  or  nonuniform  walled) ,  the  temperature 
distribution  in  the  preform,  and  the  physical  and  thermal  properties  of  materials. 

The  computer  program  begins  with  the  calculation  of  the  as-nosed  shape 
from  the  input  dimensions  of  the  nose  portion.  If  the  shell  under  consideration 
requires  a  preform  with  nonuniform  wall,  the  next  step  is  to  design  the  preform 
shape  and  divide  its  length  in  a  specified  number  of  segments.  The  actual 
simulation  begins  after  this  step.  The  die  is  moved  on  the  topmost  segment  with 
the  specified  velocity.  This  causes  the  element  to  deform  plastically.  At  this 
stage,  the  new  average  thickness  of  the  element  is  calculated.  Next,  the  average 
strains  and  strain  rate  and  the  adiabatic  heating  in  the  element  are  calculated. 
After  considering  heat  transfer  approximately,  the  new  average  temperature  of  the 
element  is  estimated.  At  these  calculated  values  of  the  strain,  strain  rate  and 
temperature,  the  flow  stress  is  determined  by  interpolation  or  extrapolation, 
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Fig  20  Functional  Flow  Chart  of  the  Program  NOSTRS 


from  a  table.  Once  the  material  flow  stress  is  obtained,  the  stress  distribution 
in  the  element  is  determined  from  the  stress  analysis.  The  meridianal  stress  at 
the  current  nose  base  is  used  to  check  for  local  bulging  and  also  to  check  for 
Euler's  buckling.  In  addition,  the  load  applied  by  the  nosing  dies  is  estimated. 

If  no  material  failure  is  predicted,  the  die  is  moved  to  cover  the  next  segment 
of  the  preform  and  all  the  above  steps  of  simulation  are  repeated  until  (a)  nosing 
is  completed,  or  (b)  material  failure  is  predicted.  In  either  case,  appropriate 
messages  are  printed  and  the  simulation  is  stopped. 

Parametric  Study 

In  order  to  illustrate  the  application  of  the  stress  and  material 
failure  analysis  and  the  associated  computer  program,  predictions  were  made  for 
cold  nosing  of  105-mm  Ml  shell,  and  hot  nosing  of  155  mm  M107  shell.  The  dimensions 
of  the  as-nosed  shape  for  155-mm  M107  and  105-mm  Ml  shell  are  taken  as  given  in 
Tables  1  and  4,  respectively. 

The  flow  stress  data  for  AISI  1045  steel  under  cold  and  hot-working 
conditions  were  taken  from  the  literature  (Ref. 25).  For  cold  nosing  of  105-mm  Ml 
shell,  the  preform  was  considered  to  have  uniform  walls,  since  the  wall  thickness 
is  not  uniform  after  nosing.  The  computer  simulation  was  carried  out  for  various 
values  of  the  coefficient  of  friction  in  the  practical  range  with  phosphate  and 
soap  lubrication.  The  load-stroke  curves  for  this  shell,  shown  in  Figure  21, 
reflect  work-hardening  of  the  shell  material  as  it  penetrates  the  die.  It  may  be 
of  interest  to  note  that  the  computer  simulation  also  predicts  considerable  warming 
of  the  shell  (150  C)  due  to  nosing. 

Since  155-mm  M107  shell  has  uniform  wall  after  nosing,  a  preform  shape 
for  this  shell  was  designed,  as  given  in  Table  3.  The  preforms  for  nosing  of 
this  shell  are  normally  preheated  by  induction  heating.  In  the  present  study,  an 
experimentally  measured  temperature  profile,  ranging  from  350  C  at  the  nose  base 
to  900  C  at  the  nose  tip  (prior  to  nosing)  was  used.  The  computer  simulation  was 
carried  out  for  various  values  of  the  coefficient  of  friction  in  the  practical 
range  with  graphite  and  water  lubrication.  The  load-stroke  curves  for  this  shell 
are  given  in  Figure  22.  The  flow  stress  of  the  material  was  taken  as  a  function 
of  local  strain,  strain  rate  and  temperature,  which  was  estimated  after  adiabatic 
heating  and  conduction  to  the  die.  As  the  shell  penetrated  the  die,  cooling  was 
greater  than  the  heating  due  to  plastic  deformation,  and  at  the  end  of  the 
simulation,  nose  tip  was  considerably  cooler  than  in  the  beginning  of  the 
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ad-Stroke  Diagram  for  Cold  Nosing  of  105-mm  Ml  Shell 
at  Various  Values  of  Friction  Coefficients 
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simulation.  This  prediction  is  in  agreement  with  observation  in  actual  hot 
nosing  of  155-mm  M107  shell.  If  the  friction  at  the  interface  was  increased 
local  bulging  was  predicted.  This  has  also  been  observed  in  practical  nosing 
operations . 


DEVELOPMENT  OF  AN  OPTIMIZATION  COMPUTER  PROGRAM 


A  comprehensive  computer  program,  named  NOSING,  was  developed  based 
on  the  analyses  of  nosing  operation  under  previous  four  sections.  Basically,  the 
computer  program  NOSING  is  an  integrated  form  of  the  computer  programs  INHEAT, 
NOSFLW  and  NOSTRS,  and  it  is  capable  of  determining  the  optimum  process  variables, 
both  in  cold  and  hot  nosing  of  shells.  For  a  given  shell,  this  program  predicts: 

(a)  As-nosed  and  preform  shape 

(b)  The  temperature  profile  in  the  preform  due  to  induction 
heating  prior  to  nosing 

(c)  The  metal  flow  during  the  nosing  operation 

(d)  The  load-stroke  curve  for  the  nosing  operation 
considering  the  flow  stress  of  the  shell  material 

as  a  function  of  strain,  strain  rate  and  temperature 

(e)  The  possibility  of  local  buckling  or  Euler's  buckling 
during  nosing. 

In  addition,  the  program  also  predicts  whether  the  nosing  can  be 
accomplished  in  single  or  multiple  hits.  A  functional  flow  chart  of  the  computer 
program  NOSING  is  shown  in  Figure  23.  As  seen  in  the  flow  chart,  the  program  is 
provided  with  options  to  bypass  intermediate  steps,  if  the  necessary  information 
is  read-in  instead  of  being  calculated. 

The  computer  program  NOSING  can  be  used  in  both  batch  and  interactive 
modes  and  requires  approximately  41,000„  words  of  memory  space.  The  input  to  the 

O 

program  is  made  through  NAMELIST  IDATA,  and  the'  input  variable  can  be  read-in 
either  in  SI  units  or  the  conventional  (inch-lb-C)  units,  depending  upon  the  user's 
option.  The  output  from  the  program  is  printed  in  units  in  which  the  input  is 
read-in.  To  illustrate  this  feature  and  to  describe  the  output  from  the  computer 
program  NOSING,  two  typical  brief  outputs  are  given  in  Figures  24  and  25.  The 
output  for  cold  nosing  of  105-mm  Ml  shell  is  in  SI  units,  as  shown  in  Figure  24, 
whereas  the  output  for  hot  nosing  of  155-mm  M107  shell  is  in  the  conventional 


Functional  Flow  Chart  of  the  Optimization 
Computer  Program  NOSING 
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units,  as  shown  in  Figure  25.  If  a  user  desires  more  detailed  output,  such  as  • 

a  printout  at  every  stage  of  simulation,  one  can  do  by  simply  specifying  an 
option  for  higher  order  of  output.  A  detailed  description  of  the  computer  i> 

program  NOSING,  including  its  various  subroutines  and  subprogram,  and  a  list 
of  important  variables  is  given  in  Appendix  E.  The  instructions  for  preparing 
input  to  the  program  are  also  included  in  Appendix  E.  By  changing  the  input 
variables  selectively,  the  output  from  the  computer  program  can  be  used  in 
determining  an  optimum  combination  of  process  variables. 


EXPERIMENTAL  EVALUATION 


In  order  to  evaluate  the  accuracy  and  capability  of  the  computerized 
mathematical  models  for  predicting  metal  flow  and  stresses  during  nosing  of  shells, 
laboratory  experiments  were  conducted.  These  cold  nosing  experiments  were  nearly 
equivalent  to  nosing  of  105  mm  Ml  shells. 

Experimental  Details 

Nine-inch  long  specimens  were  cut  from  cold-drawn  mild  steel  (AISI  1018) 
seamless  tubing  with  4-1/4-inch  0D  and  7/16-inch  thick  wall.  The  specimens  were 
annealed  and  shot  blasted  to  remove  scale.  The  remainder  of  the  scale,  developed 
during  annealing,  was  removed  by  machining  the  inner  and  the  outer  surfaces  of  the 
tubular  specimens.  The  sharp  corners  on  the  specimens  were  also  removed  by  chamfering 
both  ends.  All  the  specimens  were  marked  and  the  dimensions  (outside  diameter, 
wall  thickness  and  length)  were  measured  and  recorded.  The  cleaned  specimens 
were  then  pickled,  phosphated  and  coated  with  commercial  soap  lubricant  at  a 
cold  forging  plant  (Metal  Forge,  Columbus,  Ohio). 

The  nosing  die  was  machined  from  a  10-inch  diameter  x  11-inch  long 
piece  of  annealed  H-12  tool  steel.  The  shape  of  the  die  cavity  was  machined  to 
correspond  to  the  nosed  portion  of  a  105  mm  Ml  shell.  The  die  was  so  designed  that 
it  can  be  bolted  down  to  the  base  plate  of  Battelle's  700-ton  HPM  hydraulic  press. 

The  inner  surface  of  the  die  was  cleaned  and  polished  with  emery  paper.  The 
die  block  was  then  heat  treated  and  tempered  to  R^  52-54.  The  scale  developed 
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during  heat  treating  on  the  inner  surface  of  the  die  block  was  cleaned  by  a  wire 
brush,  and  then  the  inner  surface  was  polished  once  again  with  emery  paper. 

The  nosing  experiments  were  conducted  on  Battelle's  700-ton  HPM 
hydraulic  press.  This  press  was  instrumented  for  recording  the  load  and  the  ram 
displacement  individually  on  a  visicorder  (light-beam  oscillograph).  The  tool 
assembly  for  conducting  the  nosing  experiments  is  shown  in  Figure  26.  The  nosing 
die  was  kept  stationary  and  the  workpiece  was  pushed  in^as  shown  in  Figure  26. 

At  the  end  of  the  nosing  operation,  the  nosed  piece  was  taken  out  of  the  die  with 
the  aid  of  the  ejector. 

The  inner  surface  of  the  die  was  coated  with  MoS^  spray  prior  to  nosing. 
The  specimens  were  pushed  to  various  depths  into  the  die.  The  ram  speed  was  also 
increased  from  very  slow  (manual)  to  80  inches/minute.  The  nosing  load  and  the 
ram  displacement  were  recorded.  Several  of  the  specimens  showed  slight  buckling 
at  the  nose  base  near  the  end  of  the  stroke.  In  fact,  it  was  noticed  that  by 
increasing  the  die  penetration  by  1/8-inch  beyond  the  critical  depth  of  penetra¬ 
tion,  local  buckling  at  the  nose  base  was  initiated.  However,  since  the  shells 
elongate  due  to  nosing,  the  accurate  nose  profile  could  be  obtained  successfully 
by  depending  on  the  elongation  in  length.  At  the  end  of  nosing,  the  specimens 
were  inspected  and  the  diameter  at  the  nose  tip  and  the  final  length  of  the  specimen 
were  measured.  The  test  results  are  summarized  in  Table  5.  Figure  27  shows  a 
lubricated  specimen  prior  to  nosing  and  several  specimens  after  nosing.  Specimen 
Number  2  and  7  were  cut  along  a  diametral  plane  and  the  cross  sections  were  marked 
with  lines  at  1-inch  apart  along  the  axis,  as  shown  in  Figure  28.  The  outside 
diameter  and  thickness  of  the  nosed  specimens  along  these  lines  were  then  measured. 

Results  and  Discussion 

The  flow  stress  of  annealed  mild  steel  (AISI  1018)  was  determined  earlier 
at  Battelle  by  conducting  uniform  compression  tests.  As  shown  in  Figure  29,  the 
actual  true  stress  vs.  true  strain  curve  (o-e)  was  approximated  by  the  expression: 


<7  =  37.0  (1  +  50  e) 


0.264 


where  a  is  the  flow  stress  at  an  effective  strain  e.  The  relationship  given  by 
Equation(14)  was  used  in  the  computer  program  NOSING  to  calculate  metal  flow  and 
forces  for  the  conditions  of  laboratory  nosing  experiments. 


Bolster  Plate 


Fig  26  Tool  Assembly  for  Cold  Nosing  Experiments 


Table  5  Summary  of  Cold  Nosing  Experiments 


Specimen  No.  2  Specimen  No.  7 

Fig  28  Cross  Section  of  Specimens  2  and  7  After  Nosing 


True  Stress,  k s i 
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Fie  29  True  Stress  -  True  Strain  Curve 

for  Annealed  Mild  Steel  (AISI  1018) 


MMMHI 


71 


The  metal  flow  model  utilizes  a  constant  friction  factor  m  (0  <_  m  _<  1), 
where  the  friction  shear  stress  r  at  the  die-workpiece  interface  is  given  by 
m  a/  /3.  Since  the  value  of  the  friction  factor  at  the  die  workpiece  interface 
in  nosing  experiments  were  not  known,  predictions  were  made  for  several  values 
of  friction  factors.  As  seen  in  Figure  30,  experimentally  obtained  elongation 
in  length  fits  well  with  the  predicted  values  for  m  =  0.05.  Therefore,  for  the 
specimen  number  2,  the  predictions  of  metal  flow  were  made  using  this  value  of 
friction  factor.  As  seen  in  Figure  31,  the  predicted  change  in  wall  thickness 
and  length  of  the  shell  are  in  good  agreement  with  the  measured  values. 

The  stress  analysis,  developed  by  Nadai  (Ref.l),  utilizes  a  constant 
coefficient  of  friction  p.  Therefore,  the  load-stroke  curves  for  nosing  of 
105  mm  Ml  shell  were  generated  at  various  values  of  the  coefficient  of  friction 
using  the  computer  program  NOSING.  As  seen  in  Figure  32,  the  experimentally 
measured  curves  agree  well  with  theoretically  predicted  curve  for  p  =  0.1.  It 
is  of  interest  to  note  that  p  =  0.1  is  typically  used  in  cold  forging  analysis 
of  steel  specimens  with  phosphate  coating  and  soap  lubrication.  Further,  using 
p  *  0.1,  the  computer  program  NOSING  predicted  local  buckling  near  the  very  end 
of  the  stroke.  Slight  local  buckling  was  also  observed  in  several  specimens. 
Thus,  the  computer  program  NOSING  is  capable  of  predicting  the  load-stroke  curve 
in  cold  nosing  with  good  accuracy  and  it  is  also  capable  of  predicting  local 
buckling  due  to  nosing.  Similar  evaluations  of  the  computer  program  NOSING  for 
hot  nosing  conditions  should  be  conducted  in  the  future. 


Die  Penetration,  inch 


Fig  30  Theoretically  Predicted  and  Experimentally  Measured 
Values  of  Elongation  in  Cold  Nosing  Equivalent  to 
That  of  105  mm  Ml  Shell 


Fig  31  Theoretically  Predicted  and  Experimentally  Measured 
Nose  Profiles  (All  Dimensions  are  in  Millimeters) 


75 


REFERENCES 


1  Nadai,  A.,  "Plastic  State  of  Stress  in  Curved  Shells:  The  Forces  Required 
for  Forging  of  the  Nose  of  High-Explosive  Shells",  Forging  of  Steel  Shells, 
presented  at  the  annual  meeting  of  the  ASME,  N.Y. ,  November  29  through 
December  3,  1943. 

2  Onat,  E.  T. ,  and  Prager,  W. ,  "Nosing  of  Shells",  Technical  Report  DA798/15, 
pp  1-8,  Providence,  R.  I.,  Brown  University,  1954. 

3  Singh,  M. ,  "A  Linearized  Theory  of  Nosing  of  Shells",  Trans.  ASME,  J.  Appl. 
Mech.,  1964,  Vol.  31,  pp  535-539. 

4  Carlson,  R.  K. ,  "An  Experimental  Investigation  of  the  Nosing  of  Shells", 

Forging  of  Steel  Shells,  presented  at  the  annual  meeting  of  the  ASME,  N.  Y., 
November  29  through  December  3,  1943. 

5  Cruden,  A.  K. ,  and  Thompson,  J.  F. ,  "The  End  Closer  of  Backward  Extruded  Cans", 
NEL  Report  No.  511,  National  Engineering  Lab.,  Glasgow,  1972. 

6  Veth,  R.  D.,  Boulger,  F.  W. ,  Hull,  H.  E.,  Goble,  P.  H. ,  and  Shaw,  H.  L., 

"An  Engineering  Survey  and  Analysis  Pertaining  to  Methods  of  Manufacturing 
Artillery  Shell",  Final  Report  to  Frankford  Arsenal,  Contract  No.  DA33-019- 
507-ORD-8,  May,  1959. 

7  Pipes,  L.  A.,  "Applied  Mathematics  for  Engineers  and  Physicists",  McGraw-Hill, 
1958. 

8  Reddick,  H.  W. ,  and  Miller,  F.  H. ,  "Advanced  Mathematics  for  Engineers", 

John  Wiley,  1955. 

9  Warren,  A.  G.,  "Mathematics  Applied  to  Electrical  Engineering",  Chapman  and 
Hall,  1958. 

10  McLachlan,  N.  W.,  "Bessel  Functions  for  Engineers",  Oxford  Press,  1955. 

11  Stansel,  N.  R. ,  "Induction  Heating",  McGraw-Hill,  1949. 

12  Simpson,  P.  G. ,  "Induction  Heating",  McGraw-Hill,  1960. 

13  Rodigin,  N.  M. ,  "Induction  Heating  of  Steel  Articles  with  Normal-Frequency 

Currents",  State.  Scientific  Technical  Publishing  House  of  Literature  on 
Ferrous  and  Nonferrous  Metallurgy,  Moscow,  1950  (in  Russian). 

14  Tudbury,  C.  A.,  "Basics  of  Induction  Heating",  Forging  Industries  Associations, 
Cleveland,  Ohio,  1960. 

15  Vaughan,  J.  T. ,  and  Williamson,  J.  W. ,  "Design  of  Induction-Heating  Coils  for 
Cylindrical  Non-Magnetic  Loads",  AIEE  Transactions,  Vol.  64,  August,  1945, 

pp  587-92. 


76 


16  Vaughan,  J.  T.,  and  Williamson,  J.  W. ,  "Design  of  Induction-Heating  Coils 
for  Cylindrical  Magnetic  Loads",  AIEE  Transactions,  Vol.  65,  1946,  pp  887-92. 

17  Baker,  R.  M. ,  "Design  and  Calculation  of  Induction  Keating  Coils",  AIEE 
Transactions,  Vol.  76,  Part  II,  March  1957,  pp  31-40. 

18  Baker,  R.  M. ,  "Induction  Heating  of  Moving  Magnetic  Strip",  AIEE  Transactions, 
Vol.  64,  1945,  pp  184-189. 

19  Kothmann,  R.  E. ,  "Induction  Heating  of  a  Ferromagnetic  Cylinder  through  the 
Curie  Temperature",  Research  Report  68-1D8-IHEAT-R1 ,  Westinghouse  Research 
Laboratories,  Pittsburgh,  Pennsylvania,  April,  1968. 

20  Personal  Communication  with  Dr.  R.  E.  Kothman,  Westinghouse  Research 
Laboratories,  Pittsburgh,  Pennsylvania. 

21  Kolodii,  B.  I.,  "Temperature  Distribution  and  Thermal  Stresses  in  an 
Inductively-Heated  Hollow  Cylinder",  Prikladnaya  Mekhanika,  Vol.  5,  No.  10, 
1969,  pp  35-41. 

22  Podstrigach,  Ya.  S. ,  Goryacheva,  Z.  I.,  Kolodii,  B.  I.,  Plyatsko,  G.  V.,  and 
Komarov,  A.  N. ,  "Temperature  Fields  and  Thermal  Stresses  Produced  by  Unsteady 
Induction  Heating  of  Hollow  Cylinders",  Fiziko-Khimicheskaya  Mekhanika 
Materialov,  Vol.  6,  No.  2,  1970,  pp  113-115. 

23  Hill,  R. ,  "A  General  Method  of  Analysis  for  Metalworking  Processes",  J. 

Mech.  Phys.  Solids,  Vol.  11,  1963,  p  305. 

24  Lahoti,  G.  D. ,  and  Kobavashi,  S.,  "On  Hill's  General  Method  of  Analysis 
for  Metalworking  Processes",  Int.  J.  Mech.  Sci.,  Vol.  16,  1974,  pp  521-540. 

25  Altan,  T. ,  and  Boulger,  F.  W. ,  "Flow  Stress  of  Metals  and  Its  Application 
in  Metal  Forming  Analysis",  Trans.  ASME,  J.  Engrg.  Industry,  November  1973, 
pp  1009-1019. 


APPENDIX  A 


HP-67  PROGRAMS  FOR  AS-NOSED  AND  PREFORM  SHAPE  CALCULATION 


Two  programs,  one  for  calculation  of  as-nosed  shape  and  another  for 
designing  of  preform  shape  prior  to  nosing,  were  developed  on  an  HP-67 
programmable  hand  calculator.  These  programs  can  be  stored  on  magnetic  strips 
and  can  be  reloaded  easily. 

The  first  program  calculates  the  inner  and  the  outer  profile  of  as-nosed 
shape  from  the  nose  dimensions.  The  program  is  capable  of  treating  the  following 
two  cases: 

(a)  When  the  nose  profile  of  a  shell  is  described  by  the 
outside  ogive  radius  (R^)  and  the  uniform  wall  thick¬ 
ness  (hQ).  In  such  a  case,  it  calculates  outside  and 
inside  radii  of  the  nosed  portion  at  axial  distances 
from  the  nose  base  in  increment  of  1  inch.  First  the 
axial  distance  is  flashed,  then  the  outside  and  the 
inside  radii  of  the  nosed  shape  are  flashed  on  the 
calculator  display. 

(b)  When  the  nose  profile  of  a  shell  is  described  by  the 
outside  ogive  radius  (R^)  and  the  inside  ogive  radius 
(R^) .  In  such  a  case,  the  outside  surface  is  calculated 
as  before.  The  inside  surface  is  defined  by  another  set 
of  axial  distances  and  corresponding  inside  radius. 

First,  outside  surface  coordinates  (axial  distance  and 
outside  radius)  and  then  the  inside  surface  coordinates 
(axial  distance  and  inside  radius)  are  flashed  on  tne 
calculator  display. 

The  listing  of  this  program,  the  user's  instructions,  ar.d  the  program 
description  are  given  at  the  r nd  of  this  Appendix. 

The  second  program  in  this  Appendix  calculates  the  preform  shape  in 
order  to  obtain  a  certain  shape  after  nosing.  The  calculations  start  at  the 
base  of  the  nose  and  the  preform  wall  thickness  is  calculated  at  axial  distances 
In  increment  of  1  inch.  The  axial  distance  from  the  base  of  the  nose  is  flashed 
first  on  the  display  of  the  calculator  for  five  seconds,  and  then  the  preform 
wall  thickness  at  that  location  is  calculated  and  also  flashed  for  five  seconds. 


A- 2 


t. 


r 


At  the  end  of  the  preform  design  calculations,  if  the  coefficient  of  friction 
(y)  and  the  average  flow  stress  (a)  are  stored  in  registers  5  and  6,  respectively, 
and  the  R/S  key  of  the  calculator  is  punched,  the  program  also  calculates  an 
estimated  nosing  load. 

The  list  of  this  program,  the  user's  instructions,  and  the  program 
description  are  also  given  at  the  end  of  this  Appendix. 

Preform  Design 

Let  us  consider  an  element  of  length  dx  at  a  distance  x  of  the  preform, 
which  becomes  a  segment  of  length  ds  after  nosing.  If  the  shell  does  not  elongate 
due  to  nosing,  the  constancy  of  volume  of  the  element  yields  (Figure  A-l) : 

2ir  r  h  dx  =  2ir  r  h  ds  ,  (A-l) 

oo 

where  h  and  h  are  the  initial  and  final  wall  thicknesses,  and  the  radius  r  is 
o 

given  by: 

r  =  r  -  R,  (cosa  -  cosa) 
o  1  o 

Since  ds/da  =  cos(a  -  a  ),  from  Equation  (A-l) ,  we  have: 


h  =  — — - ,  -i  (A- 2) 

o  r  cosa(a-a  ) 

o  v  o 

If  the  shell  elongates  considerably,  an  element  originally  at  x  goes  to 
x  +  E  after  the  nosing,  where  E  is  the  axial  displacement  of  the  element.  There¬ 
fore,  r,  h,  and  ds  must  be  measured  at  (x  +  E)  when  r  and  h  are  still  measured 

o  o 

at  x.  Since 

ds  =  d(x  +  E)  cos(a  -  a  )  , 

o 

where  (a  -  a^)  is  the  slope  of  the  nose  profile  at  (x  E)i  we  have: 


k 


BGse 


(A- 3) 


A- 4 


h  =  - - -  (  1  +  ii  ) 

o  r  cos (a  -  a  )  dx 

o  o 


If  we  assume  that  5  varies  as  a  power  function  of  x  given  as: 


?  "  emax 

where  e  is  the  maximum  extension  occurring  at  the  tip  of  the  nose,  then  we  have: 
max 


hr 


h  -  ,  , 

o  r  cos(a  -  a  ) 
o  o 


r  i  ,  max 

[1  +  - —  X  n-1  ] 


(A-4) 


The  constants  n  and  e  in  the  above  expression  must  be  estimated  from  experience. 

max 

At  least  one  of  them  can  be  estimated  rather  closely  as  follows: 

Since  at  the  tip  of  the  shell,  the  state  of  pure  compression  exists  in 
the  circumferential  direction,  both  the  axial  and  the  radial  strains  are  equal  to 
half  of  the  circumf erential  strain,  £q.  The  maximum  circumferential  strain  at 
the  tip  of  the  nose  is  given  from  reduction  in  diameter  and  twice  as  much  as  the 
axial  strain  as  below: 


r  -i  r 
1  o 


d£ 


_  e  =  - -  =  - - -  (  1  +  ^  )  -  1 

r  r  cos(a.  -  a  )  dx  x=c, 

u  o  1  o  1 

~  n  e 

2  (  1  +--SSS  ) 


cos(a.  -  a  ) 
1  o 


*1 


max 


On  N  _  1 

2  '  '-u£>wl  J0/  ‘ 


Thus,  if  e  is  known,  h  is  given  bv  Equations  (A-4)  and  (A-5) 
max  o 


(A-5) 
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User  Instructions 

PROGRAM  NOSE 

CALCULATION  OF  AS-NOSED  SHAPE 

B1  <V  R2  B2 


INSTRUCTIONS 


Load  Side  1  &  Side  2  of  Card 


Input  the  known  values 


•  Radius  of  outside  ogive  curvature 


•  Distance  of  center  of  R^  from 


Nose  base  (see  Figure  1) 


•  Radius  of  inside  ogive  curvature 


•  Distance  of  center  of  R?  from 


nose  base  (see  Figure  2) 


•  Outside  radius  at  nose  base 


•  Length  of  nosed  portion 


•  Wall  thickness  at  nose  base 


or  uniform  wall  thickness  in  nose 


Calculate  outside  radius  and 


inside  radius  at  various 


axial  locations  along  the  nose 


Note:  If  R?  and  B9  not  available,  insert 


zeros  for  them.  Then  program  calculates 


inner  and  outer  surface  separatel 


Output  Form:  If  R?  ^  0,  B?  ^  0 


INPUT 

DATA/UNITS 


OUTP 
OATA/U 


etc  _ 


Output  Form:  If  R- 


x 


0  -  0. 


2 

_  _ 

_ 

3 

-  - 

— 

etc 

z  c 


Program,  Listing 


hcjij 


Program  Description 


Program  Title  CALCULATION  OF  NOSED  SHAPES 
Name  Goverdhan  D.  Lahoti 

Address  Battelle's  Columbus  Labs,  505  King  Avenue 
Citv  Columbus  ; 


Date  Aug  1,  1977 


Zip  Code  43201 


Program  Description,  Equations,  Variables,  etc. 


r  =  r  -  R  (cosa  -  cos 
o  1  o 


r  =  (r  -  h  )  -  R„  (cosa 
l  o  o  2 


l 

°  '  t/J  I 


a  =  sin  ^  (B,  /R, ) 
o  1  i 


6,  -l 

- 1  a.  =  sin  (B0/R„) 

-w.  t  1 

.  -i  ,  Bi  +  *  . 
“i  S1"  (  > 


Fig  A-2 


-1  B2  +  x 
°*2  =  sin  (  ) 


Operating  Limits  and  Warning*  When  R^ ,  E^,  »  «*ua  are  specified,  the  program  computes 

x,  r  and  r^,  in  that  order.  When  and  are  specified  0,  (when  not  available),  the 
program  computes  x,  r  and  x,  r.  pairs  for  outer  and  inner  profile. 


r*'*-  *Jk‘  .•>» 

.  MSf  iws  &&&£■ 


User  Inst  ruel  ions 


PROGRAM  PREFORM 

NOSING  OF  SHELLS  -  PREFORM  DESIGN  AND  LOAD 


INSTRUCTIONS 

INPUT 

DATA/UNITS 

Load  side  1  and  2  of  card 

Input  the  known  values 

Radius  of  “ogive  curvature 

Outside  radius  of  shell 

r0(inch) 

Length  of  nosed  portion 

2.  (inch) 

Distance  of  center  of  ogive  arc  from  the  base 


of  the  nose  (see  Figure  1)  |b(inch) 


Estimated  elongation  due  to  nosin 


Calculate  initial  wall  thicknesses  at  various 


distances  from  the  base  of  the  nose 


After  design  of  preform,  calculation  stops. 


Input  values  for  load  calculation 


Coefficient  of  friction  u 


Average  flow  stress  of  material  o(psi) 


Calculate  nosing  load 


iL-TJ  L _ I 

[  i  r  i 
[zznzz] 
nzjczz] 
r:::j  lt.:j 

i::;j  i:~j 


OUTPUT 
DAT  4,  UNITS 


KEY  ENTRY 


f  LBL  A 


KEY  CODE 


31  25  11 


A-10 

COMMENTS 


f  LB 


STO  3 


h  RTN 


f  LB 


STO  4 


h  RTN 


f  LBL  E 


STO  5 


h  RTN 


LBLf  a 


STO  6 


h  RTN 


LBLf  b 


RCL  5 


RCL  1 


33 

04 

35 

22 

31 

25 

15 

33 

05 

35 

22 

32 

25 

11 

33 

06 

35 

22 

32 

25 

12 

34 

05 

34 

01 

81  1 

Input  R 

Input  r 

c 

Input  Z 
Input  h 

Input  b 

Input  e 


(b+iO/R 

3, 


RCL  3 


h  RCI 


GTO  f  e 


DSP  4 


f  -x 


RCL 


32 

25 

13 

34 

03 

35 

34 

32 

81 

22 

31 

15 

23 

04 

31 

84 

(b+x)/R 

x 


050  |  RCL  7 


f  COS 


+ 


STEP  KEY  ENTRY  KEY  CODE 


CL  4 


3 

l 

4 

h 

5 

'd(u) 

<D 

3 

P 

✓ - V 

Q  1 

S3 

S4 

S5 

|S6 

AD-A082  106 


UNCLASSIFIED 


BATTEULE  COLUMBUS  LABS  OH  F/G  19/1 

DEVELOPMENT  of  a  COMPUTERIZED  mathematical  model  for  the  hot/co— ETC CU) 
SEP  78  6  D  LAHOTI t  T  L  SUBRAMANIAN»  T  ALTAN  DAAA25-76-C-0*27 


ARSCD-CR-78019 


Program  Listing 


STEP  KEY  ENTRY  KEY  COOE 


COMMENTS 


A- 11 

STEP  KEY  ENTRY  KEY  COOE 


32  64|  6 


COMMENTS 


(“l  -  “o> 


cos  (a  +  20) 
o 


■MRP 


mmmmmrn 


mmm 


wmm 


Program  Description 


Program  Title  Design  of  Preform  and  Calculation  of  Load  in  Nosing  of  Shells 

Name  Dr.  G.  D.  Lahoti,  Research  Scientist  Date  05/02/77 

Address  Battelle's  Columbus  Laboratories,  505  King  Avenue 


Columbus 


State  Ohio 


Zip  Code  43201 


Program  Description,  Equations,  Variables,  etc. 


(1)  When  e  =  0,  no  elongation 
max 


o  r  cos  (a  -  a  ) 
o  o 


■■  \  \ 

A  \ 


I  \ 

"T -  _  <KI 


Fig  A- 3 

l  1  ri~r, 

If  e  is  known,  n  ■  — —  (1  -  —  •  - - 

max  *  e  2  r 

max  o 

...  ^  _  hr  ..  nemax  n-1. 

or,  n  *  - ; - r-  (l  + -  x  ) 

o  rQCOS  (a  -  a  )  ^n  ' 

(3)  The  nosing  Load  P  is  given  as 


where  r  =*  r  -  R  (cos  a  -  cos  a) 

!  O  O 

f 

,  *  sin  (b/R) 

,  -1  ,  b+x. 

a  .  Sin  (  — ) 

y~4“(2)  When  e  f  0,  with  elongation 
L  max 

„ 

0  ho  •  r  -co/fa  -  .  )  «  +  S  > 

o  o 

where  e  »  e  (£)n 
max  2. 

r 

1  o,  ,  . 

— - — )  cos  (a  -  a  )  -  1 

ro  0 


P  -  2trRha  [cos(ao  +  28)  -  ey(al  "  V  cos^  +  28)  } 


where  8  *  tan-  u 


Operating  Limits  and  Warnings  (2)  When  length  of  the  nosed  portion  is  given  in  a  full  number 

and  a  fraction,  the  preform  calculations  are  carried  out  up  to  the  full  number  only,  and 

the  thickness  at  the  tip  should  be  determined  by  extrapolation.  (2)  When  e  t  0  the 

max 

initial  thickness  at  the  tip  of  the  nose  cannot  be  calculated.  When  (l  «  e  )  is  a  full 

max 

number,  the  calculated  value  of  hQ  at  the  tip  of  the  nose  should  be  ignored.  Correct 
value  can  be  obtained  by  extrapolation. 
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TEMPERATURE  DISTRIBUTION  IN  INDUCTIVELY  HEATED  TUBULAR  COMPONENTS 


APPENDIX  B 


TEMPERATURE  DISTRIBUTION  IN  INDUCTIVELY  HEATED  TUBULAR  COMPONENTS 


Heat  Generation  Due  to  Induction  Heating 

Magnetic  Field  Intensity 

In  induction  heating,  the  field  intensity  created  by  the  alternating 
current  induces  the  heating  current  in  the  charge.  Therefore,  the  distribution 
of  field  intensity  gives  a  direct  indication  of  the  distribution  of  heat  and 
temperature  gradient  in  the  charge. 

Neglecting  the  minor  variations  due  to  end  effects,  flux  variations 
and  flux  leakages,  the  magnetic  flux  intensity  surrounding  the  charge 
satisfies  Maxwell's  equation 

d2H  ,  1  dH  iHwu4ff  „ 

Jl  +  7d? - p - °  »  (B-l) 

where  H  ■  magnetic  field  intensity  surrounding  the  workpiece,  oersteds 
r  *  radius  of  test  element.  Figure  B-l,  cm  (b  £  r  <_  a) 
a  =  outer  radius  of  the  charge,  cm 
b  ■  inner  radius  of  the  charge,  cm 

u>  ■  angular  velocity  of  the  alternating  current,  radians/sec 
0  •  resistivity  of  the  material  charged,  ohm-cm 
U  *  permeability  of  the  material  charged,  gauss/oersted 

The  solution  of  Equation  (B-l)  gives  the  value  of  H  at  any  point  inside  the 
charge.  Its  general  solution  is 

H  «  A  I  (krvT  )  +  B  Ko(kn/T  )  ,  (B-2) 


*  Numbers  in  parentheses  designate  the  references  at  the  end  of  the  appendix. 


B-2 


2 

,  ,_2  4ituu  8  it  fii 

where  k  -  - -  -  - — 

P  P 

I  ,K  *  Bessel's  functions  of  first  and  second  kind  and  of  order  zero 
o  o 

A,B  ■  constants  to  be  evaluated  from  boundary  conditions. 

At  the  outer  surface,  where  r  ■  a,  flux  intensity  is  maximum  and  is  equal  to 

H  .  Therefore, 
o 


Hq  -  A  IQ(a)  +  B  KQ(a)  .  (B-3) 

For  simplicity  in  writing,  I^krv'T  )  and  similar  expressions  are  abbreviated 
as  IQ(r),  etc.,  in  Equation  (B-3)  and  subsequent  equations.  At  the  inner 
surface,  where  r  ■  b,  no  electro-motive  force  (emf)  is  induced. 

_  _  |  __P  _3H 

Er*b  ”  IrP|r-b  “  "  4tt  3r 


r«b 


(B-4) 


where  I  ■  current  density  along  the  test  element  at  a  radius  r. 
From  Equations  (B-2) ,  and  (B-4) . 


3H 

3r 


r-b 


-  kvT  ^  A  I^(b)  +  B  K^(b)J  »  0 


(B-5) 


Solving  Equations  (B-3)  and  (B-5)  simultaneously  yields: 


A  =*  -  H  K  (b)/[K  (a)  I  (b)  -  I  (a)  K  (b)] 
0  0  0  0  o  o 

B  -  H  l’(b)/[K  (a)  I '(b)  -  I  (a)  K'(b)] 

0  0  0  0  o  o 


(B-6) 


* 


Power  Loss  in  Heating 

The  heat  generated  in  the  test  element  at  a  radius  r,  thickness  3r 
and  length  31,  Figure  B-l,  can  be  calculated-  by  considering  the  total  current 
in  the  element.  The  total  current  in  the  element  is  given  by: 


x  or  o i 

r 


1  3H  i  ro 

-  7—  —  Sr  5 H 
4ir  3r 


-  ^  WT  [A  I^(r)  +  B  K^(r) ]  Sr  61 


The  resistance  of  the  test  element  is 


par 


Sr  SI 


where  a  «  angle  subtended  by  the  element  at  the  center. 


The  power  loss  or  the  heat  generation  rate  in  the  element  is  given  by: 

2 

p  -  I*R  ergs /sec 
2 

•  k  7  par  [ I A  I  (r)  +  B  K  (r) | 2 ]  Sr  SI  . 
16*2  0  0 


Substituting  (B-6)  in  (B-9) 


2  2  f 
H  k  par  I  (r)  - 
o  o _ 

16tt2  Va>  " 


*  2 

♦  K0(r) 

$  K  Ta)  6r  6<l  er«s/sec  * 


where  <p  ■  I  (b)/K  (b) 
o  0 

I  (r)  -  I  (kr»T  )  -  b«r  kr  +  i  bei  kr 
o  o 

K  (r)  ■  K  (kr/i  )  ■  ker  kr  +  i  kei  kr 
o  o 

l'(r)  -  I ' (krvT  )  -  (ber'  kr  +  i  bei'  kr)/*T 
o  o 

K  (r)  ■  K  (krvT  )  ■  (ker'  kr  +  i  kei'  kr )//T 
o  0 

H  -  0.4ir  N  I  /I 
o  c  c  w 

I  ■  rms  current  in  the  induction  coil,  amp 
c 


N  *  number  of  turns  in  the  coil 
c 


i  ■  length  of  the  coil,  cm 


B-5 


Denoting  Q  = 


I  (r)  -  *  K  (r) 
o  o 


I  (a)  -  d>  K  (a) 
o  o 


and  substituting  the  values  for 


and  k  in  Equation  (B-10)  yields: 


P  = 


4it  N  I 
c  c 

1 l 

v 


-9 


pfctr  Q  5r  Si  10  watts. 


(B-lOa) 


If  the  charge  is  heated  for  a  duration  of  5t  seconds,  energy  expended  in 
heating  the  element  is  P$t  watt  seconds,  or  joules. 


Temperature  Rise 


Thus,  the  temperature  rise,  d9,  in  the  element  during  the  time  5t 
is  obtained  as: 


d6  =  , 

Cdrad  rd l 


(B-ll) 


where  C  =  heat  capacity  of  the  material  charged,  joule/kg-K 

3 

d  »  density  of  the  material  charged,  kg/m  . 


Substituting  (B-lOa)  into  (B-ll) 


d9 


4ir  N  I 
c  c 


Hf  qtt  io-3c 


(B-lla) 


As  one  can  readily  observe,  d6  is  a  function  of  r  only.  That  is,  temperature 
gradients  are  in  the  radial  direction  only. 


B-6 

Heat  Transfer  by  Conduction,  Convection  and  Radiation 

In  order  to  predict  the  local  temperatures,  it  is  necessary  to 
consider  the  heat  equalization  due  to  heat  transfer.  Because  of  the  com¬ 
plexity  of  the  phenomenon,  it  is  appropriate  to  use  the  finite  difference 
method  for  solving  the  heat  transfer  equations.  Thus,  it  will  be  possible 
to  take  into  account  the  heat  loss  due  to  conduction  to  the  unheated  por¬ 
tion  of  the  tube,  and  due  to  convection  and  radiation  to  the  environment. 


Generalized  Difference  Equations  for  Heat  Transfer 

For  the  purpose  of  heat  transfer  analysis  by  the  finite  difference 
method,  a  portion  of  the  vertical  cross  section  of  the  tube,  shown  in  Figure 
B-2,  is  assumed  to  be  divided  into  a  number  of  trapezoidal  elements.  The 
outside  surface  of  the  tube  assumed  to  be  straight,  and  it  may  or  may  not 
be  tapered.  The  inner  surface  may  be  either  straight  or  curved  defined  by  « 

a  set  of  X,  Y  coordinates  of  points  on  the  curve.  When  the  inner  surface  is 
curved,  intermediate  points  needed  for  grid  system  are  obtained  by  polynomial 
interpolation. 

Given  these  geometric  entities,  the  length  of  the  tube  being  heated,  the 
thickness  at  the  nosing  end,  HI  in  Figure  B-2,  the  thickness  at  the  reference  loca¬ 
tion,  HO,  the  taper  angle,  4,  and  the  number  of  divisions  in  the  radial  and  axial 
directions,  a  grid  system  of  trapezoidal  elements  is  generated.  The  curved 
portion  of  the  inner  surface  of  the  tube  in  each  element  is  approximated  to  with  a 
straight  line,  as  shown  in  Figure  B-2. 

With  reference  to  Figure  B-3,  during  a  time  interval  At,  heat  trans¬ 
fer  due  to  conduction  takes  place  between  the  central  element  "O",  and  the  adja¬ 
cent  elements  "H,!  (left),  "u"  (up),  "r"  (right),  and  "d”  (down).  The  temperature 
change  in  the  volume  element  "0"  after  conduction  during  a  time  interval.  At,  is 
given  as  the  heat  balance  as  follows: 
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where,  in  addition  to  the  dimensions  given  Figures  B-2  and  B-3, 

9J  . ,9j  .  *  temperatures  of  the  element  (i,j)  before  and  after 

9  J  *  J 

time  interval  At,  respectively 

k  ,k.,k.,k  »  thermal  conductivities  between  the  central 

u  d  £  r 

element  "0"  and  the  surrounding  elements  "u", 

"d",  "l",  "r",  respectively 
RC,  .  =  centroid  of  the  central  element  "0" 

i>  j 


At  =  time  interval  during  which  heat  transfer  takes 
place 

p,C  =  Specific  gravity  and  specific  heat  of  the  material 
in  the  element  (i,j) 

h.j  ■  (aRi,j +  4Ri+i,j +  aEi,j-i +  aRi+i.i-i)/4 

ARi,j  =  (ARi,j  +  ARi+l,j)/2 


Simplifying  Equation  (B-12)  ,  the  temperature  Q'  .  of  the  element 


"0"  after  time  At  is  given  by: 
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where 


K  =  AR.  RC  cC. 

i,1  1,1 


Stability  Criterion 


The  coefficient  of  0,  cannot  be  negative.  This  case  would  be 

^  9  j 

physically  absurd,  since  it  would  mean  that  the  warmer  the  center  element  "O' 
is  before  At,  the  colder  it  would  be  after  At.  Hence, 
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The  above  inequality  determines  the  maximum  allowable  value  of 
At  for  given  thermal  conductivities  and  the  element  sizes.  This  is  called 
the  stability  criterion  or  the  convergence  condition,  and  must  be  satisfied 
in  order  to  eliminate  rapid  and  uncontrollable  computing  errors. 

Estimation  of  Element  Properties 
at  Various  Boundaries 


Depending  upon  the  location  of  a  volume  element,  the  boundary 
conditions  and  the  thermal  and  physical  properties  of  that  element  would 
vary.  For  example,  an  element  which  is  in  contact  with  air  would  also  be 
subject  to  convective  and  radiative  heat  transfer. 


Elements  in  the  Interior  of  Tube:  Referring  to  Equation  (B-13)  , 

we  have: 

kl  *  kr  '  ku  *  kd  "  V  K  '  4Ri,j  RC1,.1  'bS,  > 


where 


V  V  Cb 


thermal  conductivity,  specific  gravity,  and  specific 
heat  of  the  tube  material,  respectively. 


Elements  on  the  Front  End  of  the  Tube:  The  appropriate  difference 
equation  for  heat  balance  in  the  element  "0",  shown  by  the  broken  lines  in 
Figure  B-4a,  is  the  following: 
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Various  Volume  Elements  Considered  for 
Estimating  the  Heat  Transfer  Coefficients 
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where,  in  addition  to  the  quantities  defined  earlier, 

h^“  film  heat  transfer  coefficient  at  the  flat  end  surface 
of  the  tube 

h  *  coefficient  of  radiation 
r 

0q=  the  ambient  temperature, 
and  h^,  coefficient  of  radiation,  is  given  by: 
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where  t  ■  emissivity  of  the  surface 

a  m  Stefan-Boltzmann's  natural  constant 
Tq,T^  ^  *  absolute  temperatures  of  atmosphere  and  of  the  element, 
respectively . 

Thus,  comparison  of  Equation  (B-15)  with  Equation  (B-12) 


kt  -  kr  •- 


ka  ’  ‘b 


Az 


(h,  +  h  )  R  AR  cosa. 
f  r  i»j  i> j  i»j 


(ARJ .  .  T,  .) 
i,j  i,j 


K  ’  4Ri,j  RCi,j 


“  90  ”  atmospheric  temperature. 

Elements  at  the  Corners  of  the  Tube.  For  an  element  at  the  outer 
corner  shown  by  broken  lines  in  Figure  B-4b,  the  heat  balance  yields  the 
following  substitutions  in  Equation  (B-12): 
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where  h  =  film  heat  transfer  coefficient  at  the  free  cylindrical  surface  of 
c 


the  tube. 


Similar  expressions  for  the  inner  corner  of  the  tube  will  be: 
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Elements  on  the  Inner  and  Outer  Cylindrical  Surfaces  of  the  Tube 
For  an  element  at  the  outer  cylindrical  surface,  shown  by  broken  lines  in 
Figure  B-4c,  the  following  substitutions  are  provided  by  the  heat  balance: 
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Similar  expressions  for  the  elements  on  the  inner  cylindrical 
surface  will  be: 
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ANALYSIS  OF  METAL  FLOW  IN  NOSING  OF  SHELLS 


Introduction 


In  nosing  of  shells,  the  shape  of  the  outer  surface  of  the  shell  follows 
the  die  profile  while  the  inner  surface  is  not  supported.  Thus,  during  the  nosing 
operation,  the  part  of  the  shell  within  the  die  zone  can  move  either  radially 
inward,  or  in  the  direction  of  the  tangent  to  the  die  profile,  or  in  both  the 
directions  simultaneously.  Under  plastic  state,  metals  flow  in  the  direction  of 
least  resistance;  therefore,  the  metal  flow  in  cold  nosing  is  affected  by  the 
frictional  restraint  at  the  die-workpiece  interface  and  by  the  work-hardening 
characteristics  of  the  material  being  deformed.  In  hot  nosing,  the  metal  flow 
depends  upon  the  friction,  the  temperature  of  the  workpiece  and  the  speed  of  the 
die.  In  general,  wall  thickening  increases  with  increasing  friction,  and  in  fact, 
under  severe  friction  conditions,  the  shell  may  shrink  in  length  instead  of 

i 

elongating. 

In  this  Appendix,  the  metal  flow  in  nosing  of  shells  is  analyzed  using 
the  Hill's  general  method  of  analysis This  analysis  is  valid  for  straight- 
circular  tubular  preforms  with  uniform  wall  thickness,  as  in  the  case  of  nosing  of 
105  mm  and  smaller  size  shells.  For  larger  shells,  where  the  preform  wall  thickness 
is  not  uniform,  this  analysis  can  be  applied  only  approximately. 

Hill '8  General  Method 


The  first  step  of  this  method  consists  of  choosing  a  class  of  velocity 
fields  from  which  the  best  approximation  will  eventually  be  selected.  These  velocity 
fields  must  satisfy  all  kinematic  conditions.  The  associated  stress  (given  by  the 
material  constitutive  law)  is  either  determined  uniquely,  or  at  least  to  within  a 
hydrostatic  pressure,  if  the  material  is  incompressible.  The  associated  stress 


(1)  R.  Hill,  "A  General  Method  of  Analyses  for  Metalworking  Processes",  Int.  J. 
Mech.  Phys.  Solids,  11,  1963,  p  305. 
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distribution  in  the  deformation  zone  of  a  chosen  velocity  field  will  generally  not 
satisfy  all  the  statical  requirements.  A  question  that  now  arises  is  how  to  select, 
from  the  considered  class  of  velocity  fields,  that  one  which  most  nearly  satisfies 
the  statical  requirements. 

Noting  the  converse  of  the  virtual  work-rate  principle  for  a  continuum, 
the  final  selection  criterion  is  that  "the  statical  conditions  can  be  regarded  as 
closely  satisfied,  overall,  when 
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for  a  sufficiently  wide  subclass  of  virtual  orthogonalizing  motions  w^".  The 

surface  S  of  the  deforming  zone  generally  consists  of  three  distinct  parts,  S  =  Sc  + 

S  +  S_,  where  S„  is  the  material-tool  interface,  S  is  unconstrained  and  S  is  the 
Fit*  r  1 

interface  between  the  deformation  and  the  rigid  zones.  In  Equation(C-l) ,  denotes 
the  surface  traction  computed  from  the  considered  approximating  field  cr  „  ,  On 
surface  S„,  there  is  ordinarily  at  most  a  uniform  fluid  pressure,  p.  The  frictional 
constraint  over  surface  S  is  represented  by  a  constant  frictional  shear  stress  mk 

v 

(where  k  is  the  shear  yield  strength  of  the  deforming  material  and  m  is  a  constant, 
0<«<1),  n^  is  the  local  unit  outward  normal,  and  is  a  unit  tangent  vector 
opposed  to  the  direction  of  the  relative  velocity  of  slip  in  the  approximating  field. 
Although  the  method  is  applicable,  in  principle,  to  all  types  of  friction,  the 
representation  of  constant  frictional  stress  is  used  here. 

The  orthogonalizing  family  of  velocities  Wj  to  be  selected  must  be 
sufficiently  wide  and  should  be  just  extensive  enough  to  select  a  single  approximating 
velocity  field  from  the  particular  class  constructed  to  satisfy  the  kinematic  condi¬ 
tions.  If  the  class  of  velocity  fields  were  defined  by  equations  involving  an  unknown 
function  of  just  one  position  variable,  then  the  orthogonalizing  family  must  also 
involve  an  arbitrary  function  of  one  variable. 

Once  the  family  of  velocity  fields  is  chosen,  the  calculus  of  variations 
technique  is  applied  to  Equation (C-l) ,  treating  w^  as  a  variation.  Then  we  obtain  a 
system  of  equilibrium  equations  and  boundary  conditions,  suited  to  the  particular 
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approximating  class  and  uniquely  determining  its  best  member. 


Analysis 


As  shown  in  Figure  C-l,  consider  a  shell  of  outside  radius  r  and  uniform 

o 

wall  thickness  hQ  being  nosed  by  a  die  of  ogive  radius  R,  moving  at  a  speed  U.  A 
coordinate  system  (a,0,t),  where  a  is  the  meridional  direction,  6  is  the  circum¬ 
ferential  direction  and  t  is  the  thickness  direction,  is  selected.  It  is  assumed 
that  a, 9  and  t  are  the  principal  directions  also.  The  inner  profile  of  the  nosed 
portion  is  described  by  a  function  h(a),  which  indicates  the  wall  thickness  in  the 
t  direction  at  a  given  a.  If  r  is  the  outside  radius  at  a  given  a,  a  class  of 
approximating  velocity  fields  having  the  following  components  is  selected: 
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where  r  =  R  cosa  -  a,  and  the  prime  denotes  a  derivative  with  respect  to  a.  These 
velocities  satisfy  the  condition  of  incompressibility  given  by 
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The  velocity  field  given  by  Equation (C-2)  also  satisfies  the  velocity  boundary 
conditions.  At  the  tool-workpiece  interface,  we  have 
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Fig  C-l  Configuration  of  Shell  During  Nosing 
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At  the  free  inner  surface  of  the  nosed  portion 
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Thus,  the  velocity  field  given  by  Equation(C-2)  is  kinematically  admissible. 

In  order  to  determine  the  unknown  thickness  function  h(a)  in  Equation  (C-2) , 
we  select  the  following  virtual  orthogonalizing  velocity  field: 
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where  (a)  is  an  arbitrary  function  of  a  and  has  at  least  a  continuous  first  derivative. 
The  prime  denotes  a  derivative  as  usual.  Rewriting  the  selection  criterion  given  by 
Equation (.C-l)  in  terms  of  the  ( a,3,t)  coordinate  system,  we  have 
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where  S 


SQ,  are  deviatoric  stress  components,  o  is  the  meridiana'l  stress,  V  is 

y  t  a 


the  volume  of  the  deformation  zone,  is  the  area  of  cross  section  at  the  base  of 
the  nose  and  S  is  the  area  of  contact  between  the  shell  and  the  die. 
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Substituting  Equation (C-6)  into  Equation  ( C— 7 ) ,  we  get 


o 


where  T.  and  T  are  the  forces  acting  on  the  cross  section  at  at  =  a,  and  a  =  a  , 

1  o  1  o 

respectively.  Obviously,  is  identically  equal  to  zero  and  Tq  is  proportional  to 
force  exerted  by  the  dies.  In  the  above  equation,  dV  =  (R  -  t)  r  dB  dt  dot,  and 
dSc  =  2ttt  Rda.  Upon  substituting  and  integrating  with  respect  to  3,  we  get 
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where  F(a)  =  /  (S  -  S  )  (R  -  t)  dt 
o  a  c 


G(a)  =  /  (-  S  +  Sfi)  (R  -  t)  dt  . 
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Integrating  the  first  term  on  the  left-hand  side  of  the  above  equation  and  rearrang¬ 
ing  yields 
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Since  $(a)  is  totally  arbitrary,  the  Equation(C-S)  is  satisfied  when 
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together  with 
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Since  =  0,  F(a^)  =  0.  Equation (C-9)  is  a  second-order  ordinary  differential 
equation  in  h(a)  and  can  be  solved  using  the  following  two  boundary  conditions: 
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Once  the  solution  is  available,  the  first  of  Equations (C-10)  gives  the  required  nosing 

load  T  ■  cosa 
o  o 

The  components  of  the  deviatoric  stress  S  ,  S ^ ,  and  S^  are  obtained  from 
the  Levy-Mises  flow  rule  and  the  yield  criterion.  The  yield  criterion  can  be  written 
approximately  in  the  form  of 
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where  a  is  the  material  flow  stress.  If  i  ,  e.  and  i.  are  the  strain-rate  components, 

ct  y  u 

from  the  flow  rule 


S 

ot 


£ 

a 


2  - 

3  ° 


+  T  • 
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Substituting  Equations (C-12)  and  (C-13)  in  expressions  for  F(ot)  and  G(a) ,  we 
obtain: 


F(a)  =  |  5  h  (l  +  2  ^  •  p-)  (R  -  | 


G(o) 


-I5h^  |l  -  2  ^  |*-l) 


(C-14) 


Finally,  substituting  for  F(a)  and  G(ot)  in  Equation(C-9) ,  the  differential  eauation 
for  the  thickness  function  h(a)  is  given  as  follows: 


(2R  -  h)  p-  •  +  (5R  -  3h 

da 


rr’ ’  \  dh 

(r')2'  da 


(C-15)  * 


where  r  =  R  cosa  -  a 
r ’  =  -  R  sina 
r ' 1  =  -  R  cosa 

The  boundary  conditions  are: 


h(a  )  =  h 
o  o 


(from  F(a^)  =  0) 


(C-16) 


Equation (C-15)  is  a  second  order  ordinary  differential  equation  defining 
the  wall  thickness,  h,  of  the  shell  along  the  meridional  direction  during  nosing 
operation.  With  the  boundary  conditions (C-16)  ,  it  can  be  easily  solved  using  a 
numberical  technique.  However,  since  the  boundary  conditions  are  available  at 
each  end  of  the  interval,  some  type  of  iterative  method  is  required  for  numerical 
integration  of  Equation (C-15) .  For  this  purpose,  initially  two  guesses  of  h*  (a0) 
are  made  and  the  Equation '(C-15)  is  solved  by  a  fifth  order  Runge-Kutta  Method. 
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Normally,  these  solutions  will  not  satisfy  the  second  of  the  boundary  conditions(C-l6) . 

Therefore,  a  third  guess  value  of  h'(aQ)  is  determined  using  the  first  two  solutions 

by  linear  extrapolation  such  that  the  projected  error  in  h'(c^)  is  zero,  and  the 

entire  integration  procedure  is  repeated.  This  last  step  is  repeated  with  last  two 

solutions  until  the  error  in  h'(a^)  is  within  specified  error  bounds. 

Further,  the  Equation(.C-15)  is  singular  when  a  =  0,  since  r’(a  =  a  =  0)  =  0. 

o  o 

In  this  particular  case,  a  solution  can  be  obtained  by  taking  as  a  small  positive 
quantity  instead  of  zero.  The  error  caused  by  this  approximation  is  relatively  small. 


APPENDIX  D 


ANALYSIS  OF  STRESSES  IN  NOSING  OF  SHELLS 


A  schematic  diagram  of  nosing  of  shells  is  shown  in  Figure  D-l.  The 
length  of  the  nonuniform  walled  preform  is  divided  into  a  number  of  segments. 

If  the  number  of  segments  is  sufficiently  large,  the  variation  in  thickness  along 
the  length  of  a  segment  can  be  ignored,  and  the  average  thickness  of  the  segment 
can  be  taken  as  its  uniform  thickness.  Then  Onat  and  Prager's  analysis  can  be 
applied  to  each  element.  Once  the  element  is  inside  the  die,  its  new  average 
thickness  is  calculated  from  volume  constancy,  if  elongation  due  to  nosing  is 
given. 

Let  us  consider  the  analysis  of  stresses  for  Element  1.  The  analyses 

for  subsequent  elements  will  be  the  same,  except  for  boundary  conditions  which 

will  be  described  later.  Assuming  that  the  shell  does  not  elongate  during 

nosing,  the  rate  of  deformation  in  the  meridianal  direction,  e  ,  vanishes: 

m 

e  -  0  .  (D-l) 

m 

i 

Therefore,  from  incompressibility  condition,  the  rates  of  extension,  e  and  e 

c  n 

in  the  circumferential  and  the  normal  directions,  respectively,  must  have  equal 

absolute  values,  but  opposite  signs,  the  sign  of  e  corresponding  to  that  of  a  , 

c  c 

the  circumferential  stress.  Thus, 


e  *  -  e  <  0 
c  n 


(D-2) 


Further,  considering  an  elemental  section  at  radius  r  =  R  cosa  -  a  (see  Figure 
D-l)  and  the  thickness  h,  we  have: 


e 

c 


(D-3) 


where  r  and  h  are  the  instantaneous  rates  of  change  in  r  and  h  as  this  element  is 
pushed  toward  the  apex  of  the  die.  From  Equations  (D-2)  and  (D-3),  it  is  clear 
that  the  product  rh  maintains  a  constant  value  during  the  nosing  process.  Thus, 


rh  ■  r  h 
o  o 


(D-4) 
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where  r  and  h  are  the  initial  radius  and  initial  thickness  of  the  element, 
o  o 

Therefore,  the  variation  in  the  wall  thickness  of  the  nose  is  given  by: 


h 


r  h 
o  o 

R  cosa  -  a 


(D-5) 


The  surface  traction  exerted  by  the  die  on  the  nose  is  the  normal 

pressure  p  and  the  friction  shear  stress  f  =  pp,  where  p  is  the  coefficient  of 

friction.  The  meridianal  and  the  circumferential  stresses  produced  by  these 

surface  tractions  in  the  thin-walled  nose  far  exceed  p  and  f  in  absolute  value. 

The  state  of  stress  in  the  nose  wall  can,  therefore,  be  treated  as  plane  with 

the  circumferential  stress  a c  and  the  meridianal  stress  as  the  principal 

stresses.  Since  these  principal  stresses  will  be  compressive  and  since  o  can 

c 

be  expected  to  exceed  in  absolute  value,  according  to  Tresca’s  yield  criterion, 
we  have: 


ac  =  -  a  .  (D-6) 

The  equilibrium  of  forces  in  the  meridianal  and  the  circumferential  direc¬ 
tions  give: 

~  (rh  a  )  +  ha  sina  -  rpp  =  0  , 

da  m  c 

and  rh  a^  +  hac  cosa  +  rp  =  0  .  (D-7) 

Resolving  Equation  (D-7)  to  eliminate  p,  yields: 


~r~  (rh  a  )  +  prh  a  +  Rh  a  (sina  +  p  cosa)  =  0 
aa  m  me 


(D-8.1 


The  substitution  of  h  from  Equation  (D-5)  and  a from  Equation  (D-6)  in  Equation 
(D-8)  gives: 


da 


m  .  -  sina  +  p  cosa 

—  +  pa  =  a  - - - rr— 

da  m  cosa  -  a/R 


(D-9) 


CD-I 
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Integrating  for  a  ,  we  get: 

m 


<j  =  e 
m 


-  /u  da  (  /  f  (a)  e  da  +  C 


where 


r  \  -  sina  +  u  cosa 

f  (a)  =  a  - - - rr — 

cosa  -  a/R 


Since  at  a  =  a, ,  a  =  a',  the  specified  value  of  the  meridianal  stress  (which  is 
1  m 

zero  for  element  1),  we  have: 

al 

y(d,  -  a)  ,  1  ,  .  y(a  -  a,)  .  ,  /T.  ,, 

a  *  e  1  {a  -  /  f(a)  e  1  da }  (D-ll 

m 

a 

Once  a  is  calculated  for  element  1  from  the  known  boundary  condition,  a  for 
m  m 

element  2  and  subsequent  elements  can  be  calculated. 


Strains  in  Shell  Nose 


If  e  ,  e  and  e„  are  the  strains  in  the  meridianal,  the  circumferential 
men 

and  the  normal  directions  at  a  point  where  outside  radius  is  r  and  wall  thickness 
is  h,  then: 


e  “  ds/dx 
m 


o 


(D-12 


where  ds  and  dx  are  the  final  and  initial  length  of  a  segment  in  the  meridianal 
direction.  Since  volume  remains  unchanged  during  plastic  deformation,  we  have: 


r  .  h  .  ds  r1 

r  h  dx  ’ 

o  o 

or  (e  +  1)(e  +  1)  (e  +  1)  *  1  .  (D-lj 

men 


Thus,  if  e  and  e  are  known  from  Equation  (D-12),  e  can  be  determined  from 
c  n  m 

Equation  (D-13) .  If  e  ,  e c  and  En  are  taken  as  the  principal  strains,  the 
effective  strain  is  then  given  by: 


r 


i 

► 

t 

P 

fc 


if 

i " 

1  »i 

*- 
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f\/T 


{ (e_  - 


m 


e)2  +  (£„  "  £_)2  +  (e  -  e  )2} 
c  c  n  n  m 


(D-14) 


Strain  Rates  in  Shell  Nose 


If  h  and  r  are  the  wall  thickness  and  the  outside  radius  at  a  point  in 

the  shell  nose,  the  strain  rates  in  the  meridianal  (e  ),  the  circumferential  (e  ) 

m  c 

and  the  normal  direction  (e^)  can  be  calculated  from  the  velocity  field  given  in 
Appendix  C  as  follows: 


m 


-i  {-jh  +  jl: 

R  1  2.  +  ,2 

r  h  rh 


}  , 


1  ,  r'  , 

ec  ”  R  ^  2  *  ’ 

r  h 


1  r  h' 

Cn  '  R  {  “2  1 
rh 


(D-15) 


where  h'  and  r'  are  derivatives  of  h  and  r  with  respect  to  a. 

,  and  r  =  R  cosa  -  a, 


Since 


h  r 
o  o 


h’ 


R  cosa  -  a 


R  h  r  sina 
o  o _ 

(Rcosa-ar 


-  R  sina 


The  effective  strain  rate  is  then  elven  bv: 


D-6 


Temperatures  in  Shell  Nose 


During  nosing,  heat  is  generated  in  the  shell  due  to  plastic  deformation 
and  friction  at  the  tool-workpiece  interfaces.  Simultaneously,  heat  is  transported 
with  the  moving  material  and  heat  transfer  takes  place.  Some  of  the  generated  heat 
remains  in  the  product,  some  of  it  is  conducted  to  the  die,  and  some  may  even 
increase  the  temperature  of  the  material  moving  into  the  die.  Thus,  the  general 
problem  to  be  examined  is  that  of  time-dependent  heat  flow  in  an  incompressible 
moving  medium  with  heat  generation  in  the  medium. 

In  the  present  approximate  analysis,  it  is  assumed  that  heat  is  generated 
uniformly  in  the  shell  wall  due  to  plastic  deformation  and  friction,  and  a 
steady  state  of  temperatures  is  reached  after  an  initial  period  of  transients.  Since 
the  temperature  gradients  in  the  meridianal  direction  are  not  as  large  as  in  the 
radial  direction,  heat  transfer  in  the  radial  direction  alone  is  considered.  The 
temperature  increase  due  to  heat  generated  due  to  plastic  deformation  is  given  by: 

AT  *  a  e  /Cp  ,  (D-17) 


where  C  and  p  are  specific  heat  and  density  of  shell  material.  Since  friction 
factor  is  very  small  in  shell  nosing,  the  energy  due  to  friction  at  the  tool- 
workpiece  interface  can  be  neglected.  The  temperature  decrease  due  to  heat 
conduction  to  the  die  can  be  estimated  as: 


AT' 


2K  (T'  -  Tn) _ 

log  (rd/r)  •  VpC 


(D-18) 


where  K  =  Thermal  conductivity  of  die  material 
T'  »  T  +  AT,  T  =  initial  temperature 
Tq  =  Ambient  temperature 
r^  =  Die  outside  radio’s 
At  *  Time  interval 
V  *  Volume  of  the  element. 

Thus,  temperature,  T^,  of  an  element  inside  the  die  after  time  At  is  given  by: 


D-7 


(D-19) 


This  procedure  is  repeated  for  each  element  at  time  intervals  of  At. 
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DESCRIPTION  OF  THE  COMPUTER  PROGRAM  NOSING 


The  computer  program  NOSING  is  an  integrated  form  of  three  separate 
computer  programs — INHEAT,  NOSFLW,  and  NOSTRS — developed  to  determine  the 
optimum  process  variables  in  hot  and  cold  nosing  of  shells.  In  this  inte¬ 
grated  form,  NOSING  can  predict 

(a)  As-nosed  and  preform  shape 

(b)  The  temperature  profile  in  the  preform  due  to  induction 
heating  prior  to  nosing 

(c)  The  metal  flow  during  nosing 

(d)  The  load-stroke  relationship  during  nosing,  considering 
the  flow  stress  of  the  shell  material  as  a  function  of 
strain,  strain  rate,  and  temperature 

(e)  The  possibility  of  local  buckling  or  Euler's  buckling 
during  nosing. 

Thus,  with  the  use  of  NOSING,  generating  data  to  plot  performance  curves  is 
relatively  simple,  and  the  process  optimization  for  a  given  set  of  input  conditions 
can  be  easily  established  with  a  minimum  number  of  computer  runs. 

All  routines  in  NOSING  have  been  coded  as  subroutines  or  function 
routines.  Hence,  they  can  be  used  in  a  simple  one-line  structure  or  in  a 
superposed  overlay  structure  as  depicted  in  Figure  E-l.  In  its  present  overlay 
structure,  NOSING  requires  41000g  memory  locations.  With  one-line  structure, 
this  requirement  will  increase  to  53000g  memory  locations. 

The  computer  program  NOSING  can  be  run  both  in  batch  and  in  interactive 
mode.  When  run  in  interactive  mode,  the  input  data  is  tested  for  every  foreseeable 
error  and  the  user  is  given  an  opportunity  to  modify  the  data,  if  necessary.  In 
the  batch  run,  whenever  the  input  lacks  any  vital  data,  the  program  stops  after 
printing  appropriate  error  messages. 


NOSING 

AITKN 


Overlay  Structure  of  the  Computer  Programs  Nosing 
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This  Appendix  describes  the  functions  of  the  various  subprograms  in 
NOSING.  After  a  brief  description  of  the  input  procedure  with  suitable  examples, 
sample  outputs  from  NOSING  are  presented  to  illustrate  the  output  capability 
of  NOSING.  In  the  following  section,  each  subroutine  is  described  in  alpha¬ 
betical  order  by  its  purpose,  calling  sequence,  description  of  the  variables 
transferred  as  formal  parameters,  ’’ariables  transferred  through  common  blocks, 
brief  description  of  tha  routine  and  the  names  of  the  calling,  and  the  called 
subroutines . 

In  Table  E-l,  on  page  E-31,  all  the  variables  transferred  t hr..' ugh  th-' 
COMMON  blocks  are  described  in  alphabetical  order.  Table  E-2,  on  page  E-35, 
summarizes  the  list  of  intermediate  results  written  on  a  separate  file  facilitat¬ 
ing  program  debugging,  if  necessary. 

Input  to  the  Program  NOSING 

All  input  variables  necessary  for  process  optimization  are  read  via 

the  NAMELIST  called  IDATA.  Supplying  data  by  namelists  has  one  major  advantage. 

In  namelist  inputs,  only  those  variables  that  need  changing  are  assigned  values. 

For  this  reason,  those  variables  which  have  a  fixed  value  most  of  the  time  are 

* 

assigned  default  values  by  the  system. 

Input  to  the  program  NOSING  can  be  either  in  SI  units  or  in  the 
conventional  (inch-lb-C)  units.  Output  from  NOSING  will  be  in  the  units  used 
for  input  variables.  However,  all  computations  are  performed  internally  in  SI 
units . 

The  variables  included  in  the  namelist  IDATA  are  listed  below. 

Following  the  description  of  each  variable,  units,  wherever  applicable,  are 
given  in  square  brackets.  The  default  value  of  the  variables,  if  any,  is  given 
in  parenthesis. 


*  Stored  value  for  a  variable.  If  the  user  does  not  assign  a  value,  the 
stored  value  will  be  automatically  assigned  to  that  variable. 
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AL 

AM 

AMBTEM 

AMU 

BATCH 

BILTEM 

Bl 

B2 

COTIME 

DI 

DIEOD 

DO 

EFFICN 

ELNG 


FREQ 

GESS1 

GESS2 

HETIME 

HO 


Axial  length  of  the  nosed  portion  [mm;  inch] 

Friction  shear  factor  at  the  die-workpiece  interface 

(0.1) 

Ambient  temperature  [C]  (20.) 

Coefficient  of  friction  (0.1) 

Logical  variable.  To  be  set  to  FALSE  to  run  the 
program  interactively  (TRUE) 

Billet  temperature  [C]  (AMBTEM) 

Axial  distance  of  the  center  of  arc  R1  from  nose  base 
[mm;  inch] 

Axial  distance  of  the  center  of  arc  R2  from  nose  base 
[mm;  inch] 

Cooling  time.  Time  required  to  transfer  the  billet  from 
the  induction  heating  table  to  the  nosing  press  [seconds] 
Inside  diameter  of  the  rough  turned  shell  at  the  nose 
base  [mm;  inch]  (if  DO  and  HO  are  given,  DI  =  DO  -  2  x  HO) 
Outside  diameter  of  the  nosing  die  [mm;  inch] 

Outside  diameter  of  the  rough  turned  shell  [mm;  inch] 

(if  DI  and  HO  are  given,  DO  =  DI  +  2  x  HO) 

Efficiency  of  the  induction  coil  [percent]  (50.) 

Estimated  elongation  during  nosing.  Will  be  estimated 
by  metal  flow  analysis  unless  the  user  assigns  a  value. 
Present  metal  flow  analysis  is  not  valid  for  preforms 
with  nonuniform  wall  thickness.  Hence,  assign  a  value 
to  ELNG  if  the  preform  has  nonuniform  wall  thickness. 

[mm;  inch] 

Frequency  of  the  induction  heating  current  [cycles/sec] 
First  guess  of  the  slope  of  inner  profile  at  the  nose  base 
Second  guess  of  the  slope  of  inner  profile  at  the  nose 
base 

Duration  of  induction  heating  [seconds] 

Wall  thickness  of  the  preform  at  the  nose  base,  or  wall 
thickness  of  the  preform  with  uniform  wall  [mm;  inch]  (if 
DO  and  DI  are  given,  HO  =  (DO  -  DI)/2) 
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IC  RMS  current  in  the  induction  coil.  Real  variable, 

not  integer.  [Ampere] 

IPRINT  Printing  option.  Refer  Table  E-2.  (-1) 

ITER  Maximum  number  of  iterations  to  be  performed  during 

metal  flow  analysis.  (25) 

LC  Length  of  the  induction  coil.  Real  variable,  not 

integer.  [mm;  inch]  (AL  -  ELNG) 

NC  Number  of  turns  in  the  induction  coil.  Real  variable, 

not  integer. 

NDA  Number  of  axial  segments  for  stress  analysis  and 

temperature  estimation.  (20) 

NDOASN  Maximum  number  of  iterations  to  be  performed  during 
temperature  estimation.  Useful  during  program 
debugging. 

NDR  Number  of  radial  segments  for  temperature  estimation. 

(6) 

NSIZE  Shell  size  in  millimeters. 

R1  Radius  of  curvature  of  the  outside  nose  profile. 

[mm;  inch] 

R2  Radius  of  curvature  of  the  inside  nose  profile.  [mm;  inch] 

SIUNIT  Logical  variable.  To  be  set  to  FALSE  when  the  input 
data  are  in  conventional  (inch-lb-C)  units.  (TRUE) 

TEMP  Variable  dimensioned  25.  Temperatures  in  the  preform. 

To  be  supplied  only  when  an  estimation  by  finite-difference 
temperature  analysis  is  not  desired.  If  supplied,  there 
should  be  NDA  +  1  values,  beginning  from  the  base  to 
the  tip  of  the  nose. 

UNIPFM  Logical  variable.  To  be  set  to  FALSE  when  the  preform 
has  nonuniform  wall  thickness.  When  UNIPFM  is  FALSE, 
metal  flow  analysis  will  not  be  performed.  (TRUE) 

Axial  velocity  of  nosing  die.  [mm/sec;  inch/sec] 
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Input  by  Namelists 


When  entering  data  by  namelists,  the  first  column  is  always  left  blank. 
Following  a  $  sign  in  the  second  column,  the  namelist  name,  IDATA,  is  entered 
without  any  embedded  blanks.  Following  at  least  one  blank  space  after  the  name 
of  the  namelist,  various  parameters  are  entered  and  are  equated  to  their  values 
with  commas  between  each  parameter.  The  order  in  which  parameters  appear 
within  a  namelist  input  is  immaterial.  The  namelist  data  listing  is  terminated 
by  another  $  sign.  If  the  input  exceeds  one  line  length,  the  line  is  terminated 
after  a  comma  by  entering  a  line  feed  using  the  line-feed  (LF)  key  and  not 
the  carriage-return  key.  Every  continuation  line  is  also  started  with  a  blank 
space  in  the  first  column.  The  data  listing  is  finally  terminated  by  a  $  sign. 

In  the  batch  run,  every  continuation  card  should  have  the  first  column  blank. 
Since  any  error  in  spelling  the  variable  name  is  likely  to  cause  abnormal 
termination  of  the  program  execution,  it  is  advisable  to  have  a  list  of  the 
variable  names  included  in  the  namelist  while  trying  to  run  the  NOSING  program 
interactively. 

The  examples  given  below  show  alternative  ways  of  putting  the  data 
by  namelists. 

Example  1 

b$ IDATA  DO  =  6.2,  DI  =  4.9,  AL  =  11.,  LC  =  10.5,  SIUNIT  =  F, 
b  TEMP  =  375.,  425.,  525.,  575.,  600.,  650.,  NDA  =  5$ 

(b  denotes  blank  space) 

Example  2 

b$IDATA  DO  =  6.2,  HO  =  0.65,  LC  =  10.5,  R1  =  65.5,  IPRINT  =  5, 
bTEMP(l)  =  375.,  TEMP(2)  =  425.,  TEMP(3)  =  525.,  TEMP(4)  *  575., 
b  TEMP(5)  =  600.,  650.,  NDA  =  5$ 

As  seen  from  Example  1,  the  right  side  TRUE  or  FALSE  for  the  logical  variable 
may  be  abbreviated  to  T  and  F  respectively. 

Error  Recovery  During  the  Data  Input 

Any  error  in  the  data  input,  such  as  (1)  failure  to  enter  the  $  sign 
at  the  end  of  the  data  list,  (b)  misspelled  variable  name,  and  (c)  embedded 


/lUk  *.  IM  ****** 


Elk. 
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blanks  in  the  name  of  the  variable  will  normally  return  the  program  control 
to  the  computer  operating  system,  thereby  causing  an  abnormal  termination  of 
the  job.  However,  a  special  software  package  called  "Error  Recovery  Package" 
is  included  in  the  NOSING  system,  which  recovers  the  error  and  returns  the 
control  to  a  point  specified  by  the  package.  Thus,  whenever  an  error  occurs 
while  inputting  the  data  by  namelist,  the  message  ENTER  IDATA  VARIABLES  LIST 
will  be  printed.  In  interactive  run,  the  user  can  retype  the  list  correctly. 
In  the  batch  run,  after  the  data  ends,  a  card  with  the  word  STOP  in  the  first 
four  column,  should  be  included  to  short  the  program  in  case  of  errors  in  the 
input. 


Output  from  the  Program  NOSING 


Most  outputs  of  NOSING  are  written  directly  on  the  OUTPUT  file.  However, 
when  the  variable  IPRINT  is  assigned  a  positive  number,  several  intermediate 
results  are  listed  as  described  in  Table  E-2.  To  avoid  lengthy  output  during 
interactive  run,  these  intermediate  results  are  written  on  a  separate  file 
called  TAPE3.  This  special  file  TAPE  3  will  be  automatically  disposed  to  the 
line  printer  when  the  program  execution  is  normally  terminated.  But  when 
abnormal  termination  of  program  execution  occurs  for  any  reason,  the  user 
should  dispose  TAPE3  to  line  printer  to  obtain  the  intermediate  results.  During 
the  interactive  run,  type  DISPOSE,  TAPE3 ,  PR=C  or  ROUTE,  TAPE3,  DOPR,  TID=C 
when  the  system  is  in  COMMAND  mode.  In  the  batch  run,  add  either  DISPOSE, 

TAPE3,  PR=C.  or  ROUTE,  TAPE3,  DC=PR,  TID=C.  control  card  after  an  EXIT,  card 
in  the  control  cards  stream. 

The  first  output  from  the  program  NOSING  is  the  input  data  as  read-in, 
as  seen  in  Figure  E-2.  Other  sample  outputs  as  printed  by  NOSING  are  shown  in 
Figures  E-3,  E-4,  and  E-5..  Figure  E-3  shows  the  results  after  metal  flow 
analysis.  Figure  E-4  is  a  partial  temperature  distribution  after  induction 
heating  prior  to  nosing.  Figure  E-5  shows  the  as-nosed  profile,  the  preform 
geometry,  and  the  stress  and  load  distributions. 
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Fig  E-2  Sample  Input  Data  as  Printed  by  NOSING 
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Fig  E-3  Sample  Output  After  Metal  Flow  Analysis 


E-10 


RESULTS  FROM  THE  CONVERGED  SOLUTION 
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Fig  E-3  (Continued) 
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Fig  E-4  Partial  Sample  Output  for  Temperature  Distribution  in  the 
Preform  Corresponding  to  Inout  in  Figure  E-2 
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Fig  E-5  Sample  Output  After  Stress  and  Load  Calculations 
for  Input  in  Figure  E-2 
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The  basic  functions  of  the  main  program  and  the  various  subprograms 
of  the  computer  program  NOSING  are  briefly  described  below. 

Main  Program  NOSING 

This  main  program,  NOSING,  serves  essentially  as  a  coordinating 
routine.  It  initializes  various  input  parameters,  connects  the  input/output 
files  for  interactive  use,  initializes  the  error  recovery  package,  calls 
various  subprograms  for  preform  design  and  stress  analysis  and  finally  disposes 
TAPE3  to  the  line  printer. 


Function  ANTE 

Purpose:  Calculates  the  integral  of  the  meridianal  stress 

Function  Reference:  AINTGR  =  ANTE  (Al,  A 2) 

Al,  A2  =  Lower  and  upper  limits  of  integral 

Variables  Transferred 

Through  Common*!  ALF01,  AMU,  DO 

Description:  *  Calculates  the  integral  using  the  Simpson's  rule 

Called  by:  STRESS 


Function  ANIDIS 

Purpose:  From  the  given  set  of  coordinate  data  for  points  defining 

the  distribution  curve,  this  function  routine  computes 
the  multiplication  factor  to  estimate  the  actual  amount 
of  heat  generated  at  a  given  point  along  the  axis  of 
the  shell. 

Function  Reference:  FACTOR  =  AXIDIS (J, DELZ) 

J  =  Axial  coordinate  of  the  element  under  study 
DELZ  =  Axial  distance  between  any  two  successive  elements 
along  the  axis. 


*  All  variables  transferred  through  the  COMMON  blocks  are  described  in  Table  E-l. 
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Variables  Transferred 

Through  COMMON:  NN,  XX,  ZZ 

Description:  *  Any  intermediate  value  of  XX  at  the  desired  ZZ  value 

is  obtained  by  polynomial  interpolation  using  the 
library  routine  AITKN. 

Calling  Routines:  AITKN 
Called  Bv:  INDHET 


Subroutine  BERBEI 

Purpose :  For  a  given  value  of  X,  computes  the  Bessel's  functions 

ber(X),  bei(X) ,  ber’(X),  bei’(X),  ker(X),  kei(X),  ker’(X), 
and  kei 1 (X) 

Calling  Sequence:  CALL  BERBEI (X, BERI , INDEX) 

X  =  Coefficient  of  the  Bessel's  function 
BERI  =  An  array  of  eight  elements  carrying  back  the  ber, 
bei,  ber1,  bei',  ker,  kei,  ker’,  and  kei’  values 
INDEX  =  When  INDEX  is  set  to  zero,  ber',  bei',  ker',  and 
kei'  are  not  computed 

Description:  *  All  Bessel  functions  are  computed  by  polynomial  approx¬ 

imations  using  the  equations  given  in  the  Handbook  of 
Mathematical  Functions,  U.S.  Department  of  Commerce 

*  Two  sets  of  equations  are  used;  one  for  values  of  eight 
or  less  and  the  other  for  values  greater  than  eight. 

Called  By:  BESSEL 


Function  BESSEL 


Purpose :  Calculates  the  value  of  Q  in  Equation  (B-lOa)  given  in  Appendix  B, 
where 
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Q  = 


I  (r)  -  $K (r) 
o  o 


2 


IQ(a)  -  cJKQ(a) 

I  I 


4  =  1  (b)/K  (b) 
o  o 

I  (r)  =  I  (kr*T) 
o  o 

K  (r)  =  X  (kr  /I) 
o  o 

l’(r)  -  l'(kr/i) 
o  o 

K  (r)  =  k  (kr/i) 
o  o 


bar  kr  +  i  bei  kr 

ker  kr  +  i  kei  kr 
(ber1  kr  +  i  bei' 
(ker'  kr  +  i  kei' 


kr)//i 


kr)//i 


a  »  Outer  radius  of  the  tube 
b  =  Inner  radius  of  the  tube 
r  =  Radius  of  the  element  under  study 

k2  =  (8n2fy)/p 

f  -  Frequency  of  the  induction  voltage 
Vt,  o  =  Permeability  and  thermal  conductivity  of  t 
material  charged. 

Function  Reference:  0  =  BESSEL(X,A,B,XK) 

X  =  Radius  of  the  element,  r 
A  =  Outer  radius  of  the  tube,  a 
B  =  Inner  radius  of  the  tube,  b 
XK  =  k 

a,b,r,  and  k  are  already  defined. 

Calling  Routines:  BERBEI 
Called  Bv:  INDHET 
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Subroutine  BILLET 

Purpose:  Computer  the  thermal,  electrical,  and  magnetic  properties 
of  AISI  1045  steel  at  a  given  temperature 
Calling  Sequence:  CALL  BILLET (TEMP ,ATMT EM) 

TEMP  =  Temperature  of  the  element  under  study 
ATMTEM  =  Atmospheric  temperature 
Variables  Transferred 

Through  COMMON:  AKB,  CRB,  DENSTY,  HCB ,  HFB,  IC,  NC,  LC,  PERMEB,  RHO 

Description:  *  All  material  property  data  used  in  this  subroutine,  with 

the  exception  of  permeability,  were  obtained  from  the  ASME 
Handbook  on  "Metals  Properties",  McGraw  Hill,  1954.  Thermal 
conductivity,  heat  capacity,  and  electrical  resistivity  for 
AISI  1045  steel  were  not  readily  available.  Hence,  these 
properties  for  AISI  1040  steel  are  used  here  as  close 
approximations . 

*  Thermal  conductivity  and  electrical  resistivity  at  a  given 
temperature  are  estimated  by  linear  interpolation  from  the 
stored  set  of  discrete  data. 

*  Permeability  of  magnetic  materials  vary  considerably  with 

changes  in  the  temperature,  the  magnetic  flux  density,  and 

the  size  of  the  billet.  Adequate  data  to  meet  all  the  needs 

of  the  present  program  are  not  available  from  any  known 

source  of  information.  A  rough  approximation  given  bv  the 
„  ,  32400  ... 

rormula  p  =  — r: —  +  1  is  used,  wnere  H  is  the  maximum 
H0  o 

flux  density. 

Permeability  remains  constant  at  the  saturated  value  for 
temperatures  up  to  600  C.  Beyond  the  Curie  temperature, 
which  is  768  C  for  1045  steel,  permeability  becomes  unity. 
The  transition  between  600  and  768  C  is  very  sharp  and  is 
approximated  to  a  square  function  in  this  program. 


E-19 


Calling  RouCines:  AITKN 


Called  By:  HEATGB ,  HEATRN,  INDHET 


Subroutine  DATAIN 


Purpose :  Reads  the  process  variables  and  prints  the  supplied  data. 

Prints  appropriate  warning  messages  if  the  input  is  not 
complete.  During  interactive  run,  enables  the  user  to 
re-assign  values,  if  necessary. 

Calling  Sequence:  CALL  DATAIN 
Variables  Transferred 

Through  COMMON:  AL,  AM,  AKBTEM,  AMU,  BATCH,  BILTEM,  Bl,  B2,  COTIME,  DI , 

DIEOD,  DO,  EFFICN,  ELNG,  FREQ,  GESS1,  GESS2,  HETIME,  HO, 
IC,  IPRINT,  ITER,  LC,  NC,  NDA,  NDOASN,  NDR,  NSIZE,  Rl, 

R2 ,  SIUNIT,  TEMP,  UNIPFM,  VELO 
Called  By:  NOSING 


Subroutine  DERIVA 

Purpose:  This  subroutine  supplies  the  second  order  ordinary 

differential  equation  for  variation  in  wall  thickness 
during  nosing  in  the  form  of  first  and  second  derivatives. 
This  subroutine  is  given  in  the  EXTERNAL  mode. 

Calling  Sequence:  CALL  DERIVA  (X,F,D) 

X,  independent  variable 
F(l),  dependent  variable 

F(2),  derivative  of  F(l),  with  respect  to  X 

D(l) ,  derivative  of  F(l) ,  with  respect  to  X  (output) 

D(2),  derivative  of  D(l),  with  respect  to  X  (output) 

Called  Bv:  RK5 
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Function  FSTRS 

Purpose:  Calculates  the  required  parametric  values  for  the  inter¬ 

polation  routine  AITKN 

Function  Reference:  A  =  FSTRS (MS TRN , MTEMP , T , AS  TR , AC , AM , TEMP , STR , S TRRAT ) 

MSTRN:  Size  of  the  strain  array  of  the  flow  stress  data 

(input) 

MTEMP:  Size  of  the  temperature  array  of  the  flow  stress 

data  (input) 

T:  Temperatures  at  which  discrete  flow  stress  data 

are  available  in  the  program  (input) 

ASTR:  Strain  values  at  which  discrete  flow  stress  data 

are  available  in  the  program  (input) 

AC:  Constant  C  in  the  equation 

Flow  stress  =  C  (strain  rate)m  (input)  (E-l) 
AM:  Exponent  M  in  Equation  (E-l)  (input) 

TEMP:  Stock  temperature  (input) 

STR:  Strain  in  the  stock  (input) 

STRRAT:  Strain  rate  in  the  stock  (input) 

Description:  *  Calculate  the  intermediate  parametric  values  for  use  by 
the  library  interpolation  routine  AITKN. 

*  Call  AITKN  twice  to  estimate  the  values  of  C  and  m  in 
Equation  (E-l) . 

*  Estimate  the  flow  stress  substituting  the  values  of  C  and 
m  in  Equation  (E-l). 

Routines  Called:  AITKN 
Called  By:  MATERL 

Subroutine  GRIDSY 

Purpose:  Generates  a  grid  system  for  heat  transfer  calculations. 

Also  computes  a  number  of  parametric  values  used  in  heat 
transfer  calculations. 


Calling  Sequence:  CALL  GRIDSY 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  AL2 ,  ARI,  ARJ,  CTETA,  CTETAB,  DELZ,  DI,  DO,  DRI,  IPRINT, 

MEND,  NDA,  NDR,  NEND,  OUTP,  PHI,  R,  RC,  RSLT,  SCI,  T,  TSCR 
Description:  *  Using  the  supplied  data  for  AL2,  PHI,  DO,  and  DI 

(Refer  to  Figure  E-6) ,  R^  /s  are  computed. 


*  From  the  values  of  R.  .,  ARI.  .,  ARJ.  . ,  AR.  . ,  a .  ., 

i.J  i.J  i.J  i.j  i.J 
a,  , ,  T.  .,  and  RC.  .  are  computed, 
i.j  i.J  i.J 


RC,  . 
Figure 


is  the  centroid  of  the  element  (i.j)  as  shown  in 
E-6. 


AR  .  is  the  distance  between  the  center  points  of  the 
i.J 

elements  (i.j)  and  (i— 1, j ) ;  i.e.,  AR±  j  =  Ri  j  “  Ri_i  j. 

T,  ,  =  (AR.  .  +  AR,  .  .  +  AR  .  ,)/4. 
i.j  i.j  i+l. J  i+l, J-l 

*  When  IPRINT  _>  5,  values  of  R  are  listed  on  the  results 

file  RSLT  or  TAPE3 


When  IPRINT  _>  7,  values  of  the  centroids,  RC,  are  listed 
on  the  file  RSLT. 


When  IPRINT  >_  8,  values  of  ARI,  ARJ,  and  AR  are  listed 
on  the  file  RSLT. 


Note:  The  listed  values  of  ARJ  are  the  products  of  ARJ 
and  T;  i.e.,  ARJ.  (listed)  =  ARJ  .(computed)  x  T.  .. 


When  IPRINT  >  9,  cos  a.  .  and  cos  a.  ,  are  listed  on 

i.J  i.j 


the  file  RSLT. 


Calling  Routines:  AITKN 
Called  By:  INHEAT 
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The  temperature  distribution  is  written  on  the  RSLT  file, 
after  every  cycle  if  IPRINT  is  set  to  9,  and  is  written 
after  every  10  cycle  if  IPRINT  is  set  to  any  value  between 
3  and  9. 

BILLET,  HEATRN,  INDHET,  TMPRNT 
INHEAT 

Subroutine  HEATRN 

Purpose:  Computes  the  change  in  temperature  distribution  due  to  heat 
transfer.  The  same  routine  is  used  to  estimate  the  tempera¬ 
ture  changes  due  to  heat  equalization  during  induction 
heating  and  due  to  heat  balancing  during  cooling. 

Calling  Sequence:  CALL  HEATRN 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  AKB,  AMBTEM,  ARI ,  ARJ,  BILTEM,  CRB,  CTETA,  CTETAB,  DELT, 

DELZ,  DRI,  HCB,  HFB,  INCODE,  IPRINT,  MEND,  NEND,  R,  RC, 

RSLT,  T,  TSCR 

Description:  *  As  described  in  Appendix  B,  equivalent  k^,  k^,  k^,  k^_,  9^  ^ 

0J,,  9.  .  „ ,  9.  .,,,  and  K  are  calculated  for  elements  at 

i+l,J  i, J-l 
different  locations. 

*  The  temperature  change  in  any  particular  element  after  heat 
transfer  is  estimated  by  substituting  the  equivalent  k's, 

9's,  and  K  in  Equation  B-12  of  Appendix  B. 

*  If  IPRINT  =  7  or  IPRINT  >.  9,  values  of  I,  J,  k  ,  k  ,  k,  , 

VI  fl 

k  ,  K,  cos  or  . ,  9.  9!  ,  and  three  more  local  variables 

r*  i.J  i.J  i.j 

are  listed  on  the  RSLT  file. 

In  case  of  unexpected  results,  values  of  these  variables  will 
be  helpful  to  figure  out  the  cause  of  malfunction. 

*  Temperature  values  are  transferred  back  to  the  "T"  arrays 
from  the  scratch  "TSCR"  arrays. 


Subroutine  HEATGB 


Purpose:  This  monitor  routine  computes  the  minimum  time  interval. 

At,  in  order  to  satisfy  the  stability  criterion,  Equation 
B-14  in  Appendix  B.  By  issuing  calls  to  HEATRN  and  INDHET, 
this  routine  estimates  the  heat  generation  and  equalization 
during  induction  heating  and  during  heat  balancing  while 
the  billet  is  transferred  from  the  induction  heating 
equipment  to  the  press. 

Calling  Sequence:  CALL  HEATGB 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  AKM,  AMBTEM,  ARI,  ARJ,  BILTEM,  COTIME,  CRB,  CTETA,  CTETAB, 

DELT,  DELZ,  DRI,  HCB,  HETIME,  HFB,  INCODE,  IPRINT,  LC,  MEND, 
NDOASN,  NEND,  OUTP,  R,  RC,  RSLT,  T 
Description:  *  Temperatures  in  all  grid  points  are  initialized  by 

assigning  the  value  of  the  supplied  billet  temperature. 

*  The  minimum  time  interval,  At,  is  estimated  by  satisfying 
the  stability  criterion  at  the  most  adverse  forging  condition. 

*  From  the  prescribed  "heating  time"  and  the  "cooling  time", 
and  the  computed  minimum  time  interval.  At,  numbers  of 
iterations  required  to  cover  the  heating  period  and  the 
cooling  period  are  calculated. 

*  The  temperature  distribution  during  induction  heating  is 
estimated  by  repeated  calls  to  INDHET  and  HEATRN. 

The  temperature  distribution  is  written  on  the  RSLT  file 
after  every  cycle  if  IPRINT  is  set  to  9,  and  is  written 
after  every  10  cycles  if  IPRINT  is  set  to  any  value  between 
3  and  9. 

The  temperature  distribution  at  the  end  of  the  induction 
heating  period  is  listed  on  the  output  file. 

*  The  change  in  temperature  distribution  due  to  cooling 
during  the  transfer  of  the  billet  from  the  heating  equip¬ 
ment  to  the  press  is  calculated  by  repeating  the  calls  to 
HEATRN. 
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Calling  Routines :  BILLET 
Called  Bv:  HEATCB 


Subroutine  INDHET 


Purpose:  Calculates  the  amount  of  heat  generated  in  induction 
heating 

Calling  Sequence:  CALL  INDHET 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  CRB,  DELT,  DELZ,  FREQ,  IC,  IPRINT,  LC,  MEND,  NC,  NEND, 

PERMEB,  R,  T 

Description:  *  Amount  of  heat  generated  in  each  element  is  calculated  using 
Equation  B-lla  in  Appendix  B. 


*  If  IPRINT  8,  values  of  W,  permeability,  Q,  CRB,  A9,  and 
k  are  listed  in  the  RSLT  file,  where 


W 


f  4ir  N  I  \2 
c  c 


f  At  10 


-3 


CRB  =  C  d 

A9  =  Change  in  temperature 
Nc  =  Number  of  turns  in  the  coil 
Ic  =  Current  in  the  coil,  Amp 
Lc  =  Length  of  the  coil,  mm 
f  =  Frequency  of  the  induction  voltage,  cps 
At  =  Minimum  time  interval  for  heat  transfer  calculations 
C  *  Thermal  conductivity  of  the  material  charged 
d  *  density  of  the  material  charged 
*  (8ir^uf/p) 

Calling  Routine:  AXIDIS,  BESSEL,  BILLET 
Called  By:  HEATGB 


Subroutine  INHEAT 


Purpose :  Coordinates  the  functions  of  all  routines  used  for 
temperature  estimation 
Calling  Sequence:  CALL  INHEAT 

(No  formal  parameters) 

Variables  Transferred 

AL2,  DI,  DO,  HO,  LC,  MEND,  NDA,  OUTP,  Ql,  RSLT,  SCI,  T 
GRIDSY,  HEATGB 
NOSING 


Calculates  the  flow  stress  for  AISI  1045  steel  at  different 
temperatures,  strains,  and  strain  rates  from  the  stored 
discrete  flow  stress  data.  Also  supplies  the  specific  heat, 
thermal  conductivity,  young's  modulus,  poison's  ratio,  and 
density  of  AISI  1045  steel. 

CALL  MATERL  (STRAIN,  STRRAT ,  TEM,  FSTRES) 

STRAIN  =  Local  strain  in  the  material 
STRRAT  =  Local  strain  rate  in  the  material 
TEM  =  Local  temperature  of  the  material 
FSTRES  =  Estimated  flow  stress  (output) 

Variables  Transferred 

Through  COMMON:  ANU,  DENSTY,  E,  SPHEAT,  THCON 

Description:  *  If  the  temperature  is  less  than  500  C,  flow  stress  is 

calculated  using  the  formula: 

a  *  K(I)n  (E— 2) 

where  a  =  average  flow  stress 

c  =  average  strain 

K,n  =  constants 

*  If  the  temperature  is  more  than  500  C,  flow  stress  is 
calculated  using  the  formula: 


Through  COMMON: 
Calling  Routines: 

Called  By: 

Subroutine  MATERL 
Purpose: 


Calling  Sequence: 


where  a  =  average  flow  stress 

£  =  average  strain  rate 

C,m  =  constants. 


Calling  Routines: 

Called  By: 


For  this  purpose,  true  values  of  C  and  m  are  obtained  from 
the  stored  discrete  data  by  interpolation.  Function  FSTRS 
performs  the  interpolation  using  the  library  routine  AITKN. 
FSTRS 
STRESS 


Subroutine  NOSFLW 

Purpose:  In  addition  to  coordinating  the  functions  of  all  other 

routines  in  the  metal  flow  analysis  section,  this  routine 
solves  the  differential  equation  for  thickness  function. 
Calling  Sequence:  CALL  NOSFLW 

(No  formal  parameters) 

Variables  Transferred 


Through  COMMON: 

Description : 


A,  AL,  ALFA,  ALFO,  ALF1,  AM,  Bl,  B2,  DI,  DO,  GESS1,  GESS2 , 
H,  HO,  HP,  IPRINT,  ITER,  NDA,  RR,  Rl,  SIUNIT,  ZZ 
The  second  order  ordinary  differential  equation  is  solved 
numerically  by  the  5th  order  Runge  Kutta  method.  Since  one 
boundary  condition  is  available  at  each  end  of  the  interval 
an  iterative  method  is  employed.  That  is,  the  solution  is 
started  at  the  beginning  of  the  interval  by  assuming  a 
value  of  the  second  boundary  condition.  If  this  solution 
satisfies  the  boundary  condition  at  the  end  of  the  interval 
then  the  assumed  value  at  the  beginning  is  the  actual 
value,  and  the  calculated  solution  is  the  required  solution 
Otherwise,  the  second  boundary  condition  at  the  beginning 
is  adjusted  by  linear  interpolation  until  the  boundary 
condition  at  the  end  of  the  interval  is  satisfied  by  the 
last  calculated  solution. 
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Calling  Routine: 
Called  Bv: 


RK5,  VELCTY 
NOSING 


Subroutine  NOSTRS 

Purpose :  Defines  the  as-nosed  profile  from  the  supplied  data  and 

prints  the  geometry,  calls  PREFRM  to  define  the  preform 
geometry  and  STRESS  to  estimate  loads  and  stresses  during 
nosing . 

Calling  Sequence:  CALL  NOSTRS 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  AL,  ALFA,  ALFD1 ,  ALFD2 ,  AL2 ,  Bl,  B2,  NDA,  RRI ,  RRD,  R1 , 

R2,  TEMP,  X,  XI,  Z 
Calling  Routines:  PREFRM,  STRESS 
Called  Bv:  NOSING 


Subroutine  PREFRM 

Purpose:  Defines  the  preform  geometry  from  the  known  finish  shape 
Calling  Sequence:  CALL  PREFRM 

(No  formal  parameters) 

Variables  Transferred 

Through  Common:  AL,  ALFA,  ALF01,  AL2,  Bl,  DL,  DO,  ELNG,  H,  HO,  HP,  NDA, 

NSIZE,  R,  Rl,  SCI,  SIUNIT,  TEMP,  UNIPFM,  Z 
Description:  *  Nosing  preform  is  designed  using  simple  strain  consider¬ 
ations  given  by  Nadai.  The  estimated  elongation  due  to 
nosing  is  distributed  along  the  length  by  a  power  law. 
Preform  thickness  is  then  calculated  from  volume  constancy 
requirements . 

Called  By:  NOSTRS 


Subroutine  RK5 

Purpose :  This  routine  solves  a  system  of  first  order  ordinary 

differential  equations  with  given  initial  values,  by  a 
fifth-order  Runge-Kutta  method.  It  has  an  automatic 
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control  of  step  size  in  order  to  minimize  the  number  of 
steps  necessary  ot  obtain  a  prescribed  accuracy.  The 
routine  performs  at  least  seven  evaluations  of  the 
derivative  for  each  step. 

Calling  Sequence:  CALL  RK5  (X,Y , B, EPS ,N, FI , FP) 

X  =  a  variable  corresponding  to  independent  variable 
B  =  the  size  of  interval 

EPS  =  prescribed  relative  error.  Recommended  value, 
lO-8  <_  EPS  £  1CT4 

N  =  number  of  equations  to  be  integrated 
FI  =  a  real  variable.  Set  FI  £  0  when  entering  the 
routine  the  first  time 

FP  =  the  name  of  the  routine  supplied  by  the  user  to 
evaluate  the  derivatives.  DERIVA  in  this  case. 

Y  =  an  array  of  dependent  variable,  the  solution  array 
(output) 


Note:  This  routine  is  valid  for  N  £  25. 
Called  By:  NOSFLW 


Subroutine  STRESS 

Purpose:  Calculates  the  stress  and  the  load  during  nosing  operation. 

Calling  Sequence:  CALL  STRESS 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  ALFA,  ALOAD,  AMBTEM,  AMU,  ANU,  DENSTY,  DIEOD,  DL,  DO,  E,  EPS, 

EPSRT,  H,  HO,  HP,  IPRINT,  NDA,  PI,  R,  RRI ,  RRD,  R1 ,  SIGC,  SIGN, 
SIGT,  SIUNIT ,  SPHEAT,  TEMP,  THCON ,  VELO 
Description:  *  The  stress  analysis  used  here  was  given  by  Nadai  and 

modified  by  Onat  and  Prager.  Basically,  the  nosing  process 
is  simulated  in  a  finite  number  of  discrete  steps.  The 
stresses  at  each  step  are  calculated  in  each  segment, 
beginning  from  the  free  tip  of  the  nose.  At  each  step, 
strains,  strain  rate,  and  temperatures  are  also  calculated 
in  each  segment  and  the  flow  stress  of  a  segment  is  taken 
as  a  function  of  the  strain,  strain  rate,  and  temperature. 
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Also,  the  local  buckling  or  Euler's  buckling  is  checked 
at  each  step.  In  the  end,  a  load-displacement  table  is 
generated . 

Routines  Called:  ANTE,  MATERL 
Called  By:  NOSTRS 


Subtoutine  TMPRNT 


Purpose :  Lists  the  temperature  distribution  with  an  appropriate 
caption  on  a  prescribed  file. 

Calling  Sequence:  CALI  TMPRNT (K, PROCODE, FILNAM) 

K  =  Number  of  iterations  completed 

PRCODE  =  A  code  number  to  identify  the  caption  to  be 
printed 

FILNAM  =  File  number  on  which  the  results  are  to  be 
printed 


Variables  Transferred 

Through  COMMON:  MEND,  NEND,  T 

Description:  *  The  temperature  distribution  at  the  end  of  iteration  "K" 

is  listed  on  the  file  "FILNAM"  with  a  caption  designated 
by  "PRCODE". 

Called  By:  HEATGB 


Subroutine  VELCTY 


Purpose :  This  subroutine  calculates  the  velocity  field  in  the 

deformation  zone 
Calling  Sequence:  CALL  VELCTY 

(No  formal  parameters) 

Variables  Transferred 

Through  COMMON:  ALFA,  DO,  H,  HO,  HP,  NDA,  PI.  RR,  R1 ,  SIUNIT,  VELD,  ZZ 

Description:  *  Based  on  the  assumed  velocity  field,  the  meridianal, 

the  circumferential,  and  the  normal  components  of  velocity, 
strain  rate,  and  strain  are  calculated  for  an  assumed  thick¬ 
ness  profile  of  the  nosed  position. 

Called  By:  NOSFLW 
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DESCRIPTION  OF  THE  VARIABLES  INCLUDED 
IN  THE  COMMON  BLOCKS 


There  are  ten  labeled  common  blocks  in  NOSING.  In  the  description 
given  below,  variables  included  in  all  these  ten  blocks  are  pooled  together 
and  listed  in  alphabetical  order.  The  number  code  given  in  parenthesis 
following  each  variable  name  refers  to  the  common  block  in  which  the  particular 
variable  appears.  Explanation  of  the  number  code  is  given  below.  Default 
value  of  each  variable,  if  any,  is  given  in  brackets  below  the  variable  name. 
Units,  wherever  applicable,  are  also  given  after  the  description  of  each  variable. 


Code 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Common  Block  Name 
NOSNG 
PRFORM 
MTLFLW 
INDUCH 
SCRACH 
FLOW 
PFSTR 
HEAT 
BILPRO 
HETDIS 


Table  E-l  Description  of  Variables  in  Common  Blocks 


A(6) 

AL(1) 

AKB(9) 

ALFA(6) (7) 

ALF0(3) 

ALF01 (7) 

ALF02 (7) 

ALF1(3) 

AL0AD(7) 

AL2 (7) 

[AL  -  ELNG] 

AM(1) 

[0.01] 

AMBTEM(l) 

[20] 

ANU(2) 

ARI (8) 

ARJ(8) 

BATCH (1) 
[TRUE] 

BILTEM(l) 

[AMBTEM] 

Bid) 

B2  (1) 

COTIME (A) 
CRB ( 9 ) 


Radial  distance  of  ogive  center  from  shell  axis. (mm;  inch) 
Axial  length  of  the  nosed  portion. (mm;  inch) 

Thermal  conductivity  of  the  shell  material. (watt/m-K) 
Angular  location  of  a  segment. (radian) 

Angular  location  of  nose  base. (radian) 

Angular  location  of  nose  base  at  outside  diameter. (radian) 

Angular  location  of  nose  base  at  inside  diameter. (radian) 

Angular  location  at  the  tip  of  the  nose. (radian) 

Nosing  load  at  various  die  penetrations. (N,  Lb) 

Axial  length  of  the  nosed  portion  in  the  preform. (mm; 
inch) 

Friction  shear  factor  at  the  die-workpiece  interface. 

Ambient  temperature. (C) 

Poisson's  ratio  of  the  shell  material. 

Average  radius  in  I-direction.  Refer  to  Figure  E-6. 

Average  radius  in  J-direction.  Refer  to  Figure  E-6. 

Logical  variable.  Set  to  false  while  running  the  program 
in  the  interactive  mode. 

Initial  billet  temperature. (C) 

Axial  distance  of  the  center  of  arc  R1  from  the  nose  base, 
(mm;  inch) 

Axial  distance  of  the  center  of  arc  R2  from  the  nose  base, 
(mm;  inch) 

Time  required  to  transfer  the  heated  billet  to  the  nosing 
press. (seconds) 

Product  of  heat  capacity  and  density  of  the  shell  material. 

(Joule/m^ , K) 
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CTETA(8) 

CTETAB(8) 

DELT(8) 

DELZ(8) 

DENSTY (9) 
DI(1) 

[DO  -  2  x  HO] 
DIEOD(2) 

DL  (7) 

DO(l) 

[DI  =  2  x  HO] 
DRI (8) 

E(2) 

EFFICN(4) 

[50] 

ELNG(2) 

EPS (7) 

EPSTRC7) 

FL0W(7) 

FREQ (4) 

GESS1(3) 

GESS2 (3) 

H(6)  (7) 

HCB(9) 

HETIME(4) 

HFB(9) 

HO(l) 

[(DO  -  DI)/2] 
HP (6)  (7) 

IC(4) 


Table  E-l  (Continued) 


Cosine  of  angle  a.  Angle  at  is  shown  in  Figure  E-6. 

Cosine  of  angle  a.  Angle  a  is  shown  in  Figure  E-6. 

Maximum  time  interval.  At,  for  heat  transfer  calculations. 

Axial  distance  between  two  successive  elements.  Refer  to 
Figure  E-6. 

3 

Density  of  the  shell  material.  (Kg/m  ) 

Inside  diameter  of  the  rough  turned  shell  at  the  nose  base, 
(mm;  inch) 

Outside  diameter  of  the  nosing  die. (mm;  inch) 

Incremental  displacement. (mm;  inch) 

Outside  diameter  of  the  shell. (mm;  inch) 

AR  shown  in  Figure  B-3  of  Appendix  B. 

Young's  modules  of  the  shell  material. 

Efficiency  of  the  induction  coil. (percent) 

i 

Elongation  of  tube  during  nosing. (mm;  inch) 

Effective  strain. 

Effective  strain  rate. (1/sec) 

2 

Flow  stress. (N/mm  ,  psi) 

Frequency  of  the  induction  voltage. (cycles/second) 

First  guess  of  the  slope  of  inner  profile  at  the  nose  base. 

Second  guess  of  the  slope  of  inner  profile  at  the  nose  base. 

Instantaneous  wall  thickness  during  nosing. (mm,  inch) 

Film  heat  transfer  coefficient  at  the  cylindrical  surface 
of  the  tube. (watt/m^.K) 

Heating  time. (seconds) 

Film  heat  transfer  coefficient  at  the  flat  end  of  the  tube, 
(watt/m^.  K) 

Wall  thickness  of  the  preform  at  the  nose  base,  or  wall  thick 
ness  of  the  preform  with  uniform  thickness  wall. (mm;  inch) 

Wall  thickness  of  the  preform. (mm;  inch) 

Current  in  the  induction  coil. (Amperes) 
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Table  E-l  (Continued) 


IPRINT(l) 


ITER(l) 
[25  J 

LC  (4) 

MEND (8) 


NC(4) 

NDA(l) 

[20] 

NDOASN(l) 


NDR(l) 

[6] 

NEND(8) 

NN(10) 

NSIZE(l) 

OUTP(8) 

PERMEB(9) 

PHI (8) 

PI(1) 

[3.141592654] 

Q1,Q2,Q3,Q4(5) 

R<7)(8) 


RC(8) 

RHO(9) 


A  code  used  to  print  a  number  of  intermediate  results. 
Values  listed  under  the  direction  of  this  code  are 
summarized  in  Table  E-2.  When  the  end  results  are  not 
satisfactory,  the  intermediate  results  will  be  very  helpful 
to  figure  out  the  cause  of  program  malfunction. 

Maximum  number  of  interactions  to  be  performed  during  metal 
flow  analysis. 

Length  of  the  induction  coil  (mm). 

Maximum  coordinate  value  of  an  element  in  the  radial 
direction.  Refer  to  Figure  E-6. 

Number  of  turns  in  the  induction  coil. 

Number  of  axial  segments  for  stress  analysis  and  temper¬ 
ature  estimation. 

A  variable  used  to  limit  the  number  of  iterations  performed 
in  heat  transfer  calculations.  When  NDOASN  is  assigned  a 
value,  the  calculated  number  of  iterations  is  ignored  and 
the  heat  transfer  calculations  are  repeated  NDOASN  times. 
Helpful  while  debugging. 

i 

Number  of  radial  segments  for  temperature  calculations 

Maximum  coordinate  value  of  an  element  in  the  axial  direc¬ 
tion.  Refer  to  Figure  E-6. 

Size  of  the  arrays  XX  and  ZZ. 

Shell  size  in  millimeters. 

Name  of  the  output  file. 

Permeability  of  the  shell  material. (Gauss/cersted) 

Angle  tp  shown  in  Figure  E-6. 

Value  of  tt  used  in  area,  volume,  and  angle  calculations. 


Scratch  variables.  Used  to  transmit  values  between  routines 
in  different  segments. 

During  preform  calculations,  radius  of  as-nosed  shape  at 
various  axial  locations.  During  heat  transfer  calculations, 
radial  distance  of  the  center  of  the  element.  Refer  to 
Figure  E-6. 

Centroid  of  the  element. 

Thermal  resistivity  of  the  shell  material  (ohm-cm) . 
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RR(6) 

RRI (7) 
RR0(7) 
RSLT(8) 

Rl(l) 

R2(l) 

SIGC(7) 

SIGM(7) 

SIGT(7) 

SIUNIT(l) 

[TRUE] 

SPHEAT ( 2 ) 

T(8) 

TEMP (2) 

THC0N(2) 

TSCR(8) 

UNIPFM(l) 

[TRUE] 

VELD(l) 

X(7) 

XI  (7) 

XX(10) 

Z(7) 


Table  E-l  (Continued) 


Outside  radius  in  the  nosed  portion,  (mm;  inch) 

Inside  radius  in  the  nosed  portion. (mm;  inch) 

Outside  radius  in  the  nosed  portion,  (mm;  inch) 

Name  of  a  temporary  file,  TAPE3,  on  which  intermediate 
results  are  listed  when  IPRINT  is  set  to  a  value  between 
0  and  9. 

Radius  of  curvature  of  the  outside  nose  profile. (mm;  inch) 

Radius  of  curvature  of  the  inside  nose  profile,  (mm;  inch) 

Circumferential  stress  (N/mm  ;  psi) 

2 

Meridianal  stress  (N/mm  ;  psi) 

2 

Normal  pressure  (N/mm  ;  psi) 

Logical  variable.  Set  to  FALSE  when  the  input  data  are 
lb-inch-C  units. 

Specific  heat  of  the  shell  material.  (J/Kg) 

Temperature  array.  In  the  subprogram  GRIDSY,  T  is  used 
as  a  scratch  array 

An  array  of  temperature  distribution  in  the  preform  prior 
to  nosing. 

Thermal  conductivity  of  the  shell  material. (W/m.K) 

A  scratch  array  to  store  intermediate  temperature  values. 

Logical  variable.  Set  to  FALSE  when  the  preform  has  non- 
uniform  wall  thickness. 

Axial  velocity  of  the  nosing  die.  (mm/sec;  inch/sec) 

Axial  location  on  outside  surface  of  nosed  portion. (mm;  inch) 

Axial  location  on  inside  surface  of  nosed  portion,  (mm;  inch) 

X  coordinates  of  the  points  defining  the  heat  distribution 
along  the  axis  of  the  tube. 

Axial  location  of  a  segment. (mm;  inch)  In  heat  transfer 
calculations,  Z  coordinates  of  the  points  defining  the  heat 
distribution  along  the  axis  of  the  tube.  (mm;  inch) 

In  metal  flow  analysis,  ZZ  are  the  location  of  axial 
segments. (mm;  inch) 


ZZ(6) (10) 
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List  of  Variables  Written  on  the  File  TAPE3 
for  Different  Settings  of  IPRINT 

IPRINT  is  a  variable  included  in  the  namelist  IDATA.  When  IPRINT 
is  set  to  a  value  between  0  and  9,  during  program  execution,  several  inter¬ 
mediate  results  are  written  on  a  separate  file  named  TAPE3.  At  the  end  of  the 
program  execution,  the  file  is  automatically  disposed  to  the  line  printer. 
These  intermediate  results  may  be  helpful  to  find  possible  errors  if  there  are 
any.  The  list  of  variables  written  on  TAPE  3  for  different  settings  of  IPRINT 
are  summarized  in  this  table.  When  IPRINT  is  set  to  a  value,  say  5,  all  the 
variables  corresponding  to  5,  in  the  table  shown  below,  as  well  as  those 
corresponding  to  5,  will  be  written  on  TAPE  3. 

Table  E-2 


IPRINT  Setting 


Variables  Listed 


2 

3 

5 

2 

7 


8 


*  DELT  and  number  of  iterations  required  to  estimate  temp¬ 
erature  distributions  at  the  end  of  induction  heating  and 
after  cooling  from  HEATGB. 

*  Temperature  distributions  (i.e.,  ^i,j  for  i  =  1  to  MEND 
and  j  =  1  to  NEND)  after  each  ten  iterations  during 
induction  heating  and  during  cooling  from  HEATGB. 

*  Radius  of  each  grid  (i.e.,  j  for  i  =  1  to  MEND  and 
j  =  1  to  NEND)  from  GRIDSY. 

*  ALFA,  ZZ,  RR,  H,  AND  H’  after  each  iteration  from  NOSFLW. 

*  Centroid  of  each  grid  (i.e.,  RCi^j  for  i  =  1  to  MEND  and 
j  =  1  to  NEND)  from  GRIDSY. 

*  i»j  »ku>lcl»kr>kd,K,cosai  ..cosc^  .  ,9^  .,9'^.  -s  and  three  other 
variables  for  i  =  1  to  &END  and’j  =  I  to  NliND  from  HEATRN. 

*  ARI^ j ,ARit j ,ARJi> j  for  i  =  1  to  MEND  and  j  =  1  to  NEND 
from' GRIDSY. 

Note:  Values  listed  under  ARJj_?j  are  actually  the 
product  of  ARJisj  and  T^j. 
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IPRINT  Setting 


Table  E-2  (Continued) 


Variables  Listed 


*  W,u  (permeability),  Q,  C,  A6,  and  K  for  each  value  of  I  and 
J  from  INDHET  where 


W  = 


/4tt  Nil2  , 

C  C’  f  6t  10~3 


*  C  =  Thermal  capacity  x  density 


P 

*  cosot^  j ,  cosa ^  ^  for  i  =  1  to  MEND  and  j  =  1  to  NEND  from 
GRIDSY. 

*  0£  j  for  i  =  1  to  MEND  and  j  =  1  to  NEND  after  each 
iteration  during  induction  heating  and  during  cooling 
from  HEATGB 

*  i,j  ,ku,kd,ki,kr,K,cosai#j  ,j  ,cosaij  ,j  .Qij  and  three  other 
variables  for  i  =  1  to  MEND  and  j  =  1  to  NEND  from  HEATRN. 


DISTRIBUTION 


Department  of  Defense 

Director  of  Defense  Research  and 
Engineering  Office 
Attn:  Mr.  S.  Persh 
Washington,  DC  20301 

Director 

Defense  Advanced  Research  Projects  Agency 
Attn:  Dr.  E.  C.  Van  Reuth 
Dr.  C.  Lehner 
Dr.  E.  Blase 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Defense  Documentation  Center 
Attn:  TIPDR 
Cameron  Station 
Alexandria,  VA  22314 

Department  of  the  Army 

Commander 

US  Army  Material  Development  Readiness  Command 
Attn:  DRCDMD 
DRCDMR 
DRCDMD- T 

DRCD-L,  Foreign  Science  5  Tech  Div 

DRCDE-E,  Eds on  Gardner 

DRCDE-W 

DRQ-rr ,  L.  Croan 
DRCMT,  Col  N  Vinson 
DRCDE-DE,  E.  Lippi 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

Aberdeen  Proving  Ground 
Attn:  STEAP-TL,  Technical  Library 
DRXSY 
DRXSY-GA 


Copies 


Director 

US  Army  Air  Mobility  Research  §  Development  Lab 
.Ames  Research  Center 

Attn:  Mr.  Paul  Yaggy  1 

Moffet  Field,  CA  94035 

Commander 

US  Army  Air  Mobility  R§D  Labs 

Attn:  SAVDL-ST  1 

Fort  Eurstis,  VA  23604 

Commander 

US  Army  Material  Development  §  Readiness  Command 
Scientific  and  Technical  Infor  Team-Europe 


Attn.  DRXST-STL  1 

APO,  New  York  09710 

Commander 

Rock  Island  Arsenal 

Attn:  DRXIB-MT  1 

Technical  Information  Div  1 

SARRI-ER  1 

Commander 

US  Army  Harry  Diamond  Labs 

Attn:  AMXDO-TIB  1 

2800  Powder  Mill  Road 
Adelphia,  MD  20783 


Commander 

US  Army  Materials  §  Mechanics  Research  Center 


Attn:  DRXMR-X,  Dr.  E.  Wright  1 

DRXMR-PT  1 

Technical  Information  Div  1 

Wratertown,  MA  02172 

Director 

US  Army  Maintenance  Management  Center 

Attn:  DRXMD-A  1 

Lexington,  KY  40507 

Commander 

Watervliet  Arsenal 

Attn:  DRDAR- LCB-TL  1 

SARWV-PPI  1 

SARWV-QA  1 

SARJW-RD  1 

SAR1W-RDR,  Dr.  T.  Davidson  1 

Technical  Information  Div  1 

Watervliet,  New  York  12139 


A 


> • 


* 


2 


* 


Copies 


Commander 

US  Army  Research  Office 
P.O.  Box  12211 

Attn:  Ur.  George  Mayer,  Director  1 

Dr.  E.  Saibel  1 

Research  Triangle  Park,  N'C  27709 

Commander 

US  Army  Natick  Research  5  Development  Command 
Attn:  DRXRE,  Dr.  E.  Sieling  1 

Technical  Information  Div  1 

Natick,  MA  07160 

Commander 

US  Army  Foreign  Science  5  Tech  Center 

Attn:  W.  Mar ley  1 

J.  Bollendorf  1 

7th  Street,  NE,  Building  220 
Charlottesville,  VA  22901 


Commander 
Edgewood  Arsenal 
Attn:  SAREA-MT 

Aberdeen  Proving  Ground,  MD  21010 
Commander 

US  Army  Mobility  Equipment  Research  §  Development 


Command 

Attn:  STSFB-MfM,  Mr.  W.  Baer  1 

DRSf‘E-RZT  1 

Ft.  Belvoir,  VA  22060 

Commander 
Redstone  Arsenal 

Attn:  Technical  Information  Div  1 

Huntsville,  AL  35809 

Commander 

Ballistic  Missile  Defense  Systems 

Attn:  BNDSC-TS  1 

DRDMI-R  1 

DRCPM-LCE,  B.  Crosswhite  1 

P.O.  Box  1300 
Huntsville,  AL  35809 


Commander 

Rocky  Mountain  Arsenal 

Attn:  Technical  Information  Div  1 

Denver,  CO  80240 


3 


Coi 


Director 

US  Army  Production  Equipment  Agency 
Attn:  DRXPE,  J.  Callaugher 
Rock  Island,  IL  61201 

Director 

US  Array  Advanced  Materials  Concept  Agency 
Attn:  Technical  Information  Div 
2461  Eisenhower  Avenue 
Alexandria,  VA  22314 

Commander 

US  Array  Tank-Automotive  R$D  Command 
Attn:  DRDTA-RKA,  V.  Pagano 
DRDTA-KP 

DRDTA-RE,  C.  Bradley 
DRDTA-Z,  J.  Panks 
Warren,  MI  48090 

Commander 

US  Army  Tank-Automotive  Materiel  Readiness  Command 
Attn:  DRSTA-E 
Warren,  MI  48090 

Commander 

US  Army  Aviation  Systems  Command 

Attn:  DRSAV-ERE 

P.O.  Box  209 

St.  Louis,  MO  63166 

Commander 

US  Array  Troop  Support  Command 
Attn:  DRSTS-PLC 
4300  Goodfellow  Blvd 
St.  Louis,  MO  65120 

Director 

USDARCOM  Intern  Training  Center 
Attn:  DRXMC-ITC-PPE 
Red  River  Army  Depot 
Texarkana,  TX  75501 


Copies 


Commander 

US  Army  Armament  R§D  Command 

Attn:  DRDAR-CG,  MG  B.  Lewis  1 

DRDAR-TD,  Dr.  R.  Weigle  1 

DRDAR-TDR,  E.  Eichelberger  1 

DRDAR-SC,  Dr.  D.  Gyorog  1 

DRDAR-LC,  Dr.  J.  Frasier  1 

DRDAR-SCM,  J.  D.  Corrie  3 

Dr.  E.  Bloore  1 

Project  File  1 

DRDAR-SCM- P,  Mr.  I.  Betz  1 

Dr.  K.  Iyer  3 

DRDAR-SCM-P,  Mr.  F.  Lee  10 

DRDAR-LCU-M  1 

DRDAR-LCU,  Mr.  Bushey  1 

DRDAR-LC  1 

DRDAR-P'tl  1 

DRDAR-PBiM-GA,  G.  O'Brien  1 

DRDAR-SC,  Mr.  W.  Dittrich  1 

Commander 

Rep,  Alabama  Army  Ammunition  Plant  1 

Childersburg,  AL  35044 

Commander 

Rep,  Badger  Army  Ammunition  Plant  1 

Baraboo,  WI  53919 

Commander 

Rep,  Burlington  Army  Ammunition  Plant  1 

Burlington,  NJ  08016 

Commander 

Rep,  Comhusker  Army  Ammunition  Plant  1 

Grand  Island,  MB  68801 

Commander 

Rep,  Gateway  Army  Ammunition  Plant  1 

St.  Louis,  MO  63143 

Commander 

Rep,  Hays  Army  Ammunition  Plant  1 

Pittsburgh,  PA  15207 


Commander 

Holston  Army  Ammunition  Plant 
Kingsport,  TN  37662 


5 


1 


Copies 

Commander 

Indiana  Army  Ammunition  Plant 

Charlestown,  IN  47111 

1 

Commander 

Iowa  Army  Ammunition  Plant 

Burlington,  IA  52600 

1 

Commander 

Rep,  Joliet  Army  Ammunition  Plant 

Joliet,  IL  60*436 

1 

Commander 

Kansas  Army  Ammunition  Plant 

Parsons,  KS  67357 

1 

Commander 

Lake  City  Army  Ammunition  Plant 

Independence,  J.D  64050 

1 

Commander 

Lone  Star  Army  Ammunition  Plant 

Texarkana,  TX  75501 

1 

Commander 

Longhorn  Army  Ammunition  Plant 

Marshall,  TX  75671 

1 

Commander 

Louisiana  Army  Ammunition  Plant 

Shreveport,  LA  71102 

1 

Commander 

Milan  Army  Ammunition  Plant 

Milan,  TN  38358 

1 

Commander 

Newport  Army  Ammunition  Plant 

Newport,  IN  47966 

1 

Commander 

Radford  Army  Ammunition  Plant 

Radford,  VA  24141 

1 

Commander 

Rep,  Ravrenna  Army  Ammunition  Plant 
Ravenna,  OH  42266 

Commander 

Riverbank  Army  Ammunition  Plant 
Riverbank,  CA  95367 

Commander 

Scranton  Army  Ammunition  Plant 
Scranton,  PA  18501 

Commander 

Rep,  St.  Louis  Army  Ammunition  Plant 
St.  Louis,  .MO  63160 

Commander 

Rep,  Sunflower  .Army  Ammunition  Plant 
Lawrence,  KS  66044 

Commander 

Twin  Cities  Army  Ammunition  Plant 
New  Brighton,  MN  55112 

Commander 

Volunteer  Army  Ammunition  Plant 
Attn:  SARVO-T 
Chattanooga,  TN  37401 

Commander 

Anniston  Army  Depot 
Attn:  DRXAN-DM 
Anniston,  AL  36201 

Commander 

Corpus  Christi  Army  Depot 

Attn:  DRXAD-EFT 

Corpus  Christi,  TN  78419 

Commander 

Fort  Wingate  Depot  Activity 
Attn:  DRXFW-M 
Gallup,  NM  87301 


/ 


Commander 

Tobvhanna  Amy  Depot 
Attn:  DRXTO-'E-B 
Tobyhanna,  PA  1S466 

Commander 
Tooele  Amy  Depot 
Attn:  DRXTE-SEN 
DRXTE-EMD 
Tooele,  UT  84074 

Commander 

Letterkenny  Army  Depot 
Attn:  DRXLE-M 
DRXLEBM 

Clambers burg ,  PA  17201 

Commander 

Lexington- Blue  Grass  Armv  Depot 
Attn:  DRXLX-SE- 1 
Lexington,  KY  40507 

Commander 

New  Cumberland  Army  Depot 

Attn:  DRXXC-S'l 

New  Cumberland,  PA  17070 

Commander 
Pueblo  Army  Depot 
Attn:  DRXPU-IE 
DRXPU-SE 

Pueblo,  CO  81001 
Commander 

Red  River  Army  Depot 
Attn:  DRXRIK'M 
Texarkana,  TX  75501 

Commander 

Sacramento  Army  Depot 
Attn:  DRXSA-M4E-LB 
Sacramento,  CA  95815 


1 

1 


1 


1 


1 

1 


1 


1 


8 


Copies 


Commander 
Seneca  Armv  Depot 
Attn:  URXSE-SE 
Romulus,  NY  14541 

Commander 
Sharpe  Army  Depot 
Attn:  DRXSH-SO 
DRXSH-M 

Lathrop,  CA  95330 

Commander 
Sierra  Armiy  Depot 
Attn:  DRXSI-DQ 
Herlong,  CA  96113 

C.  Department  of  the  Navy 

Chief 

Bureau  of  Ships 
Department  of  the  Navy 
Attn:  Code  343 
Washington,  DC 

Chief 

Bureau  of  Aeronautics 
Department  of  the  Navy 
Attn:  Technical  Information  Div 
Washington,  DC 

Chief 

Bureau  of  Weapons 
Department  of  the  Navy 
Attn:  Technical  Information  Div 
Washington,  DC  20025 

Commander 

Naval  Air  Development  Center 
Johnsville,  Aero  Materials  Dept 
Attn:  Mr.  Forrest  Williams 
Warminister,  PA  18974 

Commander 

Naval  Air  Systems  Command 
Department  of  the  Navy 
Attn:  AIR  5203,  Mr.  R.  Schmidt 
AIR  604 

Washington,  DC  20560 


1 


1 

1 


1 


1 


1 


1 


1 

1 


9 


=*1 


Copies 


Officer  in  Charge 

US  Navy  Material  Industrial  Research  Office 

Attn:  Code  227  1 

Philadelphia,  PA  19112 

Commander 

Naval  Ships  Systems  Conmand 
Department  of  the  Navy 

Attn:  Code  03423  1 

Washington,  DC  20023 

Commander 

US  Naval  Weapons  Laboratory 

Attn:  Technical  Information  Div  1 

Dahlgren,  VA  22448 

Commander 

US  Naval  Engineering  Experimental  Station 

Attn:  WCTRL-2,  Materials  Lab  1 

Annapolis,  MD  21402 

Commander 

US  Naval  Ordnance  Laboratory 

Attn:  Code  WM  1 

Silver  Spring,  MD  20910 

D.  Department  of  the  Air  Force 


Director 

Air  Force  Materials  Laboratory 

Attn:  AFML,  Technical  Library  1 

LTE  1 

LTM  1 

LTN  1 

Wright- Patterson  AF3 
Dayton,  Oh  45455 


Director 

US  Naval  Research  Laboratory 
Attn:  Mr.  W.  S.  Pellini,  Code  6300, 


Metallurgy  Div  1 

Mr.  W.  J.  Ferguson  1 

Mr.  J.  Baker  1 

Dr.  H.  C.  Posey  1 

Dr.  F.  Rosenthall  1 

Mr.  C.  Sanday  1 

Washington,  DC  20375 


10 


Director 

Naval  Ships  Research  §  Development  Center 
Attn:  Mr.  Abner  R.  V. 'illmer 

Chief  of  Metals  Research 
D.  W.  Taylor 
Code  042,  Tech  Lib 
Bethesda,  MD  20084 

Director 

Air  Force  Armament  Laboratory 
Attn:  AFATL/DLOSL 
ADTC/DLJW 

Eglin  AFB,  FL  32542 
Director 

Air  Force  Weapons  Laboratory 
Attn:  Technical  Information  Div 
Kirtland  AFB,  Nil  87118 

Director 

Air  Force  Materials  Laboratory 
Wright'Patterson  AFB 
Attn:  AFML/LLD,  Dr.  T.M.F.  Ronald 
AFML,  Tech  Library 
Dayton,  OH  45433 

Commander 

Aeronautical  System  Div 
Wright-Patterson  Air  Force  Base 
Attn:  Technical  Information  Div 
Dayton,  OH  45433 

Commander 

.Air  Research  8  Development  Command 
Andrews  Air  Force  Base 
Attn:  RDRAA 
Washington,  DC  20025 


Director 

National  Academy  of  Science 
Attn:  Materials  Advisory  Board 
2101  Constitution  Avenue,  N.W. 
Washington,  DC  20418 


11 


1  rr<  *,r? 


F 


Copies 


V' 


T 


Director 

National  Aeronautics  5  Space  Administration 
Attn :  Code  RPM 
Federal  Building  #10 
Washington,  DC  20546 

Director 

National  Bureau  of  Standards 
Attn:  Technical  Information  Div 
Washington,  DC  20025 

F.  Private  Organizations 

Battelle  Memorial  Institute 
Attn:  Metals  §  Ceramic  Information  Center 
Dr.  T.  Leontis 
505  King  Avenue 
Columbus ,  OH  43201 

Materials  Research  Laboratory,  Inc. 

Attn:  Dr.  E.  J.  Ripling 
1  Science  Road 
Glenwood,  IL  60425 


1 


1 


1 

1 


1 


i* 


V 


* 


l 


